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ABSTRACT 


This dissertation considers two different topics. In the first part of the dissertation, 
we show that a positive proportion of the primes have the property that if any one 
of their digits in base 10, including their infinitely many leading 0 digits, is replaced 
by a different digit, then the resulting number is composite. We show that the same 
result holds for bases b € {2,3,--- ,8,9, 11,31}. 

In the second part of the dissertation, we show for an integer b > 5 that if 
a polynomial f(2) with non-negative coefficients satisfies the condition that f(b) is 
prime, there are explicit bounds D4(b) and D3(b) so that if the degree of f(a) is < Da, 
then f(x) is irreducible; and if the degree of f(x) is < D3 and f(x) is reducible, then 
f(x) is divisible by the shifted cyclotomic polynomial ®,(x2 — b). Furthermore, in 
the case that b > 26, there are explicit bounds Dg(b) and D(b) so that if the degree 
of f(x) is < Dg and f(x) is reducible, then f(x) is divisible by either ®4(a — b) or 
®3(x—b); and if the degree of f(a) is < D and f(x) is reducible, then f (2) is divisible 
by ®,4(a — b), &3(a — b), or ®3(x — b). Furthermore, we show that for each b > 5 and 
each n € {3,4,6}, the bound D,, is sharp. 
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CHAPTER 1 


INTRODUCTION 


In this dissertation, we detail a variety of number theoretic results. Broadly speak- 
ing, these results fall into either the topic of digit changing or that of polynomial 
irreducibility. The results in this introduction are therefore split into two sections 
depending on which topic they pertain to, and the proofs that follow in the disserta- 
tion are split into two chapters based on the same criterion. A third section of the 
introduction lists several results which the author helped to establish, the proofs of 


which appear in other dissertations. 


1.1 RESULTS PERTAINING TO CHANGING THE DIGITS OF A NUMBER 


In 1978, M. S. Klamkin [19] proposed the following problem: Does there exist any 
prime number such that if any digit (in base 10) is changed to any other digit, the 
resulting number is always composite? Shortly afterwards, P. Erdés [5] provided a 


proof that there are an infinite number of such primes. The first several are 


294001, 505447, 584141, 604171, 971767, 1062599, 


1282529, 1524181, 2017963, 2474431, 2690201, 3085553. 


A complete list up to 10° can be found through a link in [26]. More recently, T. Tao 
[30] established that a positive proportion of the primes satisfy the conditions in 
Klamkin’s problem, and J. Hopper and P. Pollack [17] showed that a positive propor- 
tion of the primes become composite if one modifies any single digit and appends a 


bounded number of digits at the beginning or end. To clarify, following the wording 


of J. Hopper and P. Pollack [17], we define a prime p as being digitally delicate if 
it has the property given in Klamkin’s problem, that is, if any one of its digits is 
changed, then the resulting number is composite. If mg(a) denotes the number of 
digitally delicate primes up to x and, as usual, 7(a) denotes the number of primes 


up to x, then Tao’s work implies 


ma(10") 32393 


= = 0.000071031... 
(101°) 455052511 


it seems likely the liminf above is very small. 

The work of Erdés can be generalized to any base, and Tao as well as Hopper 
and Pollack gave rather general results for arbitrary bases. For example, note that a 
prime p in base b has r = r(p,b) = |log,p| + 1 base b digits. Thus, we can write a 
prime p in base b as Dee, d;b’ where 0 < d; < b—1 for each j. Tao’s work allows 
one to extend Klamkin’s problem by including a number of leading 0’s in the prime 
p. Taking s > r, we can write 
0 ihe ae 


p=) .d;b’, where d;= 
a0 some number in {0,1,...,b5-—1} if0<j<r. 


Hence, we can view p as having > r digits by allowing leading 0’s. 


Theorem 1.1 (Tao). For every integer C > 0, a positive proportion of the primes p 
have the property that in every base b € [2,C] if any of the C([log,(p)| + 1) base b 
digits of p, including (C — 1) ([log,(p) + 1) leading 0’s, is changed to another digit, 


then the the resulting number is composite. 


As an interesting related example, the first prime in base 10 for which the number 


of leading 0’s that can be included in Klamkin’s problem is greater than the number 


of digits of the prime is the prime p = 354975121. The smallest prime that occurs by 


changing any one digit of p including leading 0’s is 
70000000000354975121. 


In this dissertation, we consider primes with leading 0’s on a larger scale. 


Definition 1.2. A widely digitally delicate prime is a prime with the property that 
if any one of its digits, including any of its infinitely many leading 0’s, is replaced by 


any different digit, then the resulting number is composite. 


While this notion makes sense in any base, we will focus in Chapter 2 on widely 
digitally delicate primes in base 10. We establish the following theorem about such 


primes. 


Theorem 1.3. A positive proportion of the primes in base 10 are widely digitally 


delicate. 


We do not know of any examples of widely digitally delicate primes in base 10. 
None of the digitally delicate primes up to 10° listed through [26] are widely digitally 
delicate. 

We prove Theorem 1.3 in Sections 2.1 and 2.2. In Section 2.1, we exhibit specific 
covering systems in a manner similar to [10] to show that for sufficiently large primes 
p satisfying a certain congruence, increasing any digit (possibly a leading 0) in the 
base 10 expansion of p results in a composite number. In Section 2.2, we first use a 
partial covering similar to the argument of Erddés in [5] to show that for p satisfying 
a second congruence, most changes arising from decreasing a digit of p result in a 
composite number. We then mimic Tao’s argument in [30] by using a sieve to show 
that many of these primes become composite when one of the few remaining nonzero 


digits is decreased, so that Theorem 1.3 follows. 


Because our argument depends on exhibiting particular covering systems which 
allow for increasing a digit in base 10, our proof does not generalize directly to other 
bases. However, the rest of our argument only depends on the existence of such a 
covering system, so to prove the statement analogous to Theorem 1.3 for a given 
base 6 it suffices to exhibit a covering system which establishes that any increase 
of a digit in base b results in a composite number. We have found such covering 
systems for bases b € {2,3,...,11}, which are listed in Appendix A. However, our 
methods become difficult to replicate as b grows. We also looked at the base b = 31 
since finding a covering is simplified in the case that b — 1 has distinct small prime 
factors. We were able to show that a positive proportion of the primes are widely 
digitally delicate in base 31, and the coverings for that argument also appear in 
Appendix A. We emphasize again though that we do not know whether a result 
similar to Theorem 1.3 holds for an arbitrary base. 

The corresponding case of b = 2 in Theorem 1.3 has received prior attention in the 
literature in a different form. The theorem in this case is connected to a number being 
both a Sierpinski and a Riesel number. A Sierpifiski number is a number n such that 
n-2*+1 is composite for all non-negative integers k, and a Riesel number is a number 
n such that n-2* — 1 is composite for all non-negative integers k. D. Ismailescu and 


P. S. Park [18] have shown that if 
n = 10439679896374780276373 (mod 66483084961588510124010691590), 


then n is both a Sierpinski and a Riesel number. 

The notion of a Sierpiiski prime relates to the notion of a widely digitally delicate 
prime in the following way: Let n € Zt. One can show, as was shown for Riesel 
numbers in Lemma 4 of [8], that if there is a finite set of primes P so that for all 
sufficiently large k, some prime p € P divides n - 2* + 1 (respectively n- 2" — 1), then 
for all k in Z, the number n+2* (respectively n—2*) is divisible by some prime p € P. 


Since the proof in [18] which obtained the congruence above relied on determining 


such a set P for each of the numbers n-2* +1 and n-2*—1 where k is a non-negative 


integer, and since the set P has elements less than any element of 
{|10439679896374780276373 + 2*| : ke Z+U {oF}, 


it follows that primes p satisfying this congruence and, hence, a positive proportion 
of primes, satisfy that p+ 2" and p — 2* are composite for every non-negative integer 


k. In particular, the prime 
10439679896374780276373 


has the property that if any of its base 2 digits, including any of its infinitely many 
leading 0’s, is replaced by a different digit, then the resulting number is composite. 
The smallest proven widely digitally delicate prime in base 2 is 2131099, which we 
obtained by looking for digitally delicate Sierpinski primes. We searched the numbers 
listed through [27] for digitally delicate primes and 2131099 was the smallest such 
prime. It is widely digitally delicate since 2131099 is digitally delicate and has been 


proven to be a Sierpinski number with covering set 
P= {3,5, 7, 13,17, 2414, 


so that the numbers 2131099 + 2* are also composite for all non-negative integers k 
since any base 2 digit change in 2131099 results in a number larger than 241. 

We have also used the coverings listed in Appendix A to find a prime which is 
widely digitally delicate in base 3. Beyond this, we have been unable to find widely 
digitally delicate primes in other bases. 

There are a number of interesting open questions related to Theorem 1.3 that can 


be asked, which we address in Section 2.3. 


1.2. RESULTS PERTAINING TO POLYNOMIAL IRREDUCIBILITY 


If d,d,_1...dido is the decimal representation of a prime, then a result of A. Cohn [22] 
asserts that 


f(t) =d,2"+d,12" 1 4+---+det+dy 


is irreducible over the integers. This result can be generalized in several directions if 
we view f(x) as a polynomial with non-negative coefficients and f(10) prime. Ques- 
tions arise in this setting: Does the irreducibility of f(a) depend on the coefficients 
being less than 10? Can we find a degree n so that if f(a) has degree less than n, 
then f(x) is irreducible? Is 10 special, or do similar results hold when we replace the 
condition that f(10) is prime with the condition that f(b) is prime for some positive 
integer b 4 10? 

Various answers to these questions can be found in the literature. The result 
of Cohn has been extended to all bases b > 2 by J. Brillhart, A. Odlyzko and M. 
Filaseta [2]. In [6], M. Filaseta extended this to base b representations of kp where 
k is a positive integer < b and p is a prime, and M. R. Murty [23] has an analog in 
function fields over finite fields. Furthermore, [2] allows the coefficients d; in Cohn’s 
theorem to satisfy 0 < d; < 167 rather than 0 < d; < 9; and later M. Filaseta [7] 
showed that the d; need only satisfy 0 < d; < 10°°d,,, and further that simply d; > 0 
suffices ifn < 31. 

More recently, work by S. Gross and M. Filaseta [9] extended this line of inves- 
tigation in other ways. They showed that if f(x) is a polynomial with non-negative 
coefficients, f(10) prime, and the degree of f(x) less than or equal to 34, then f(z) is 
reducible only in the case that f(z) is divisible by ®4(2 —10) = 2? — 20x +101, where 
®,,(x) is the nth cyclotomic polynomial. They further showed that if 34 is instead 
replaced by 36, then f(x) is reducible only if it is divisible by one of ®4(a — 10) and 
®3(x — 10). 


Building on this work in [4], M. Cole, S$. Dunn, and M. Filaseta extended this 
result in the setting where f(b) is prime, 2 < b < 20, to find similar bounds D(b), 
D,(b), and D2(b) so that if the degree of f(a) is less than or equal to D(b), then 
f(x) is irreducible, if the degree of f(x) is less than or equal to D,(b) and f(z) is 
reducible, then ®4(a — b) divides f(x), and if the degree of f(x) is less than or equal 
to D2(b) and f(x) is reducible, then f(x) is divisible by ®4(x — b) or ®3(a — b). Our 
main goal for this section of the dissertation will be to extend these results via the 


theorem below. Note that we have modified the notation used from that in [4]. 


Theorem 1.4. Fiz an integer b > 5, and for n € {3,4,6} set 


T T 
— and Ey . 
a ay a " sam aI 


Let f(x) € Za] with non-negative coefficients and with f(b) prime. If the degree of 


Dy = D,(b) = 


f(x) is < Dg, then f(x) is irreducible. Additionally, if the degree of f(x) is < D3 
and f(x) is reducible, then f(x) is divisible by ®4(a—b) and not divisible by ®3(x —b) 
or g(x —b). Furthermore, in the case that b > 26, if the degree of f(x) is < De and 
f(x) is reducible, then f(x) is divisible by either ®4(x — b) or ®3(x — b) and not by 
Oo(x — b). Lastly, in the case that b > 26, if the degree of f(x) is < D and f(x) is 
reducible, then f(x) is divisible by ®4(a — b), &3(x — b), or ®3(ax — b). 


The proof of Theorem 1.4, which appears in Chapter 3, relies heavily on the 
techniques established in both [9] and [4]. Throughout the dissertation, irreducibility 
will refer to irreducibility in Z[z]. We show in Chapter 3 that the bounds D,,(b) above 
are sharp, in that they are the largest numbers having the stated property. However, 
we do not know if D(b) is sharp. In particular, the denominator in the floor in D 


corresponds to an approximation of 


arctan ( 


asi) 


As b grows, one can allow for less precise approximations of 3 in D which admit 


the same results as those in Theorem 1.4. 


On the coefficient side of this question, it was shown in [9] that if f(x) is a 


polynomial with non-negative coefficients bounded above by 
49598666989151226098104244512918 


and f(10) is prime, then f(x) is irreducible over Z. S. Gross and M. Filaseta also 


showed that if the coefficients were instead bounded above by 
8592444743529135815769545955936773, 


then f(z) is either irreducible over Z or divisible by ®4(x — b) = x? — 20x +101. Fur- 
thermore, they showed that these bounds are sharp, i.e., they exhibited polynomials 
with maximum coefficient one more than the numbers displayed above which were 
reducible (and in the second case, not divisible by x? — 20x + 101). 

This work was generalized further in [4], which found bounds M;(b) such that if 
the coefficients of f(x) are bounded above by M,(b), and f(b) is prime for an integer 
b € [2, 20], then f(x) is irreducible over Z. They also found bounds M2(b) such that if 
the coefficients of f(x) are bounded above by Mo(b), and f(b) is prime for 3 <b < 5, 
then f(z) is either reducible or divisible by ®3(a — b). Similarly, if 6 < 6 < 20, and 
the coefficients of f(a) are bounded above by Mo2(b), then f(x) is either reducible or 
divisible by ®4(a — 6). Furthermore, they established that the upper bounds M,(b) 
are sharp for 3 < b < 20, and that the upper bounds M(b) are sharp for 4 < b < 20. 


We state here an extension of the results in [4] to all integers b > 2. 


Theorem 1.5. Let b € Z with b > 2 and let n € {3,4,6}. Set 


r= 0) = laeere| 


Let f(a): = ue a;x) with a; > 0 for each j and f(b) prime. Let 


B” = max (BY, By?) 


By = | (, 2 ' (B+ Re(Gn)) P24 — ie (1 An + By) 
and 
(mn) __ Dr, —1 Dy—2k-2(__ 2\k = 
Bye s Ik = 1 (b Se Re(n)) ( Im(¢n) ; (1 An a By), 


(1.1) 


for (A3, Aa, Ag) = (2b — 1, 2b, 2b+ 1) and (B3, By, Be) = (b? —b 1? 52 b 1). 


Also, for b > 20 let 


__ (b — 1.5221)*(b — 2.5221) _ | (-1)" 
By = I a aot (7/02) with ky => | 5 (1.2) 


where @ = arctan(0.8444/(b—0.2)), and for 2 < b < 20 define By to be values chosen 
in [22, Table 10]. Then the following hold. 

For b = 2, if eacha,; < 7, then f(x) is irreducible. For 3 <b <5, if eacha; < B®), 
then f(x) is irreducible. For b > 5, if each aj < Bo, then f(x) is irreducible. For 
3<b<5, if eacha; < BO) and f(x) is reducible, then f(x) is divisible by &3(a — b). 
For b > 5, if each a; < B®) and f(x) is reducible, then f(x) is divisible by ®4(x — b). 
For b > 69, if each a; < Bi and f(x) is reducible, then f(x) is divisible by ®3(x — b) 
or ®4(a—b). For b > 69, if each a; < By, and f(x) is reducible, then f(x) is divisible 
by &3(x — b), B4(a — b), or Be(x — Bb). 


The proof of Theorem 1.5 will be discussed further in the dissertations of J. Foster 


and J. Juillerat. 


1.3. OTHER RESULTS 


We now proceed to other results of the author. In [12], J. Foster, J. Juillerat, and 


the author use Newton polygons to establish that a certain family of polynomials is 


irreducible. The work stems from [14], in which Heim, Luca and Neuhauser study 
the functions 


exp (S90 *) SRG) 


n>1 n>0 
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where g = e with 7 in the upper complex half-plane, and g : N > N is an arith- 
metic function normalised such that g(1) = 1. Given these restraints, the polynomials 


P3(zx) can be shown to satisfy the recursive relationship 
= ES 

for all n > 1, where P§(x) = 1. Of particular interest in [14] is the case that 
g(n) = a(n) = Nand. In this setting, the roots of the polynomials P?(x) appearing 
as coefficients in the resulting g-power series can be shown to dictate the vanishing 
properties of the nth Fourier coefficients of powers of the Dedekind eta function (see, 
for example, [16], [21]). 

Using the fact that n < o(n) < n?, Heim, Luca and Neuhauser seek to gain 
information about P’(«) and their roots by considering the family of polynomials 
P(x) where g(n) =n and g(n) = n?. 


In the case that g(n) = n, the recursive formula for P9(x) gives the closed form 


ae feel ee 
Pu(a)=a)) ( k Feant 
Using a result of Schur [24] regarding the irreducibility of the generalized Laguerre 
polynomials, they show that in this case when g(n) = n the polynomial P9(x) is x 
times an irreducible polynomial. However, this result does not apply when g(n) = n?, 


in which case they show P9(x) = «P,(2) where 


had 


The polynomial P,(x) is Eisenstein when n — 1 is prime, and hence P,(2) is irre- 


ducible. Heim, Luca, and Neuhauser leave as an open question whether the P,() 
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are irreducible for general n. Using Newton polygons, J. Foster, J. Juillerat, and 
the author establish the following theorem, which will be discussed further in the 


dissertations of J. Foster and J. Juillerat. 


Theorem 1.6. The polynomials P, (2) are irreducible over Q for all positive integers 


nN. 


Using methods similar to those established by Heim and Neuhauser in [15], J. Fos- 
ter and the author have generalized Theorem 1.6 to hold when g(n) = n‘t', t € Zt. 


In this case, the polynomials P9(x) take the closed form 


oni) t) kl 


where 


S(nlkt)= SD mie mt = [2") (2 1 ote? 4 3t73 4 )') 


Mit +M,p=N 


and |a"|(-) denotes the nth coefficient operator. In the case that t = 1, S(n|k, 1) 


n+j—1 
2j-1 


specializes to the binomial coefficient ( ) appearing in P, (x). 
One can show that the formula for S(n|k,t) above simplifies to the (n — k)th 
coefficient of x in the expansion 
X G . dl (Tac) Zi) » ( Ce se 
ko+---+kp_1=k ’ ’ i=0 j=0 
where the numbers A(t,7) are the so-called Eulerian numbers. Using this formula 


and an argument involving Newton polygons similar to that in [12], J. Foster and the 


author establish the following theorem. 


Theorem 1.7. Let t € Z* and let g(n) =n'*'. Then the polynomial 


stih) t) ie 


Pea) =e 


is x times an irreducible polynomial. 


CHAPTER 2 


WIDELY DIGITALLY DELICATE PRIMES 


2.1 INCREASING A DIGIT 


We remind the reader of the following definition and theorem from the introduction. 


Definition 1.2. A widely digitally delicate prime is a prime with the property that 
if any one of its digits, including any of its infinitely many leading 0’s, is replaced by 


any different digit, then the resulting number is composite. 


Theorem 1.3. A positive proportion of the primes in base 10 are widely digitally 


delicate. 


We begin by establishing that a positive proportion of primes p satisfy that the 


numbers p+ 6-10" are composite 
V6 € {1,2,...,8,9} and VkeZtu {0}. 


This then will account for changing any leading 0 of p as well as increasing any digit of 
p that is not a leading 0. The approach here is similar to that used in [10]. We achieve 
the above by constructing finite sets of primes Ps; which arise from constructing a 
covering (defined below) for the powers of k appearing in the numbers p + 6 - 10°. 
Then p will have the property we want if we choose p from an appropriate residue 


class A modulo the product of primes in P; UP: U---U Po. 


Definition 2.1. A finite system of congruences x = a; (mod m,), 1 <i < t, is called 
a covering of the integers (or simply a covering) if each integer satisfies at least one 


congruence in the system. 
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For example, one can check that the system 


x=0 (mod 2) 
x=1 (mod 3) 
x=3 (mod 4) 
x=5 (mod 6) 


x=9 (mod 12) 


is a covering of the integers. Although this example is a covering where the moduli 
are distinct, for our purposes, we do not want to require that the moduli need be 
distinct. 

We next state and prove a simple lemma, similar to Lemma 1 in [10], which 
clarifies our use of covering systems. By way of notation, we use c(p) to refer to the 


multiplicative order of 10 modulo a prime p with gced(p, 10) = 1. 


Lemma 2.2. Let A€ Z*. Fora fixed 6 € {0,...,9}, suppose we have distinct primes 


D1,---,pt, each relatively prime to 10, satisfying the following. 
(i) There exists a covering of the integers 
k=b; (mod c(p;)), 1<i<t. 
(ii) The number A satisfies each of the congruences 
A=-6-10% (mod p), 1<i<t. 
Then, for each k € Z* U {0}, the number 
A+6-10* 


is divisible by at least one of the primes p; where 1 <i <t. 
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Proof. Suppose the conditions in the lemma hold. Let k € Z*+ U {0}. By (i), there is 
ani € {1,...,t} such that k = b; (mod c(p;)). Write k = c(pi)q + );. Since c(p;) is 


the multiplicative order of 10 modulo p;, we have 10°47 = 1 (mod p;). Hence, 


A+6-10*=A+6-10% (mod ). 


From (ii), we deduce A + 6-10" =0 (mod p,), completing the proof. 


We will determine an explicit example of primes pj, ..., p: satisfying the conditions 
in Lemma 2.2 for each 6 € {0,...,9}. We will denote this set of primes {pi,..., pz} 
by Ps. The idea is to consider primes p= A (mod p;---p;) so that when a digit of p 
is increased by 0, the resulting number is divisible by some p; with 1 <i < t (as the 
lemma implies). In order for primes p to exist with p = A (mod p,---p;), we will 
require gcd(A, p,--+p,) = 1. 

We now proceed to choose for each 6 € {1,...,9} appropriate A, primes p1,..., De, 
and corresponding b;,..., 5, so that (i) and (ii) of Lemma 2.2 hold. For 6 € {2,5, 8}, 
we take A = —2 (mod 3) and 


Py = Ps = Ps = {3}. 


In other words, for 6 € {2,5,8}, we take t = 1 and p,; = 3. As c(3) = 1, we can 
take b; = 0. Observe then that (i) and (ii) of Lemma 2.2 hold. Also, with A = —2 
(mod 3), we see that gced(A, 3) = 1. This application of Lemma 2.2 can be viewed as 
a degenerate case, where the covering of the integers in (i) is the single congruence 
k =0 (mod 1). 

For our subsequent choices, we will make use of the factorizations of 10" +1 in 


[3] to find primes p so that n = c(p). For example, one checks that 


2 = c(1l), 4 = c(101), 


8. = lis), -28. = e137). 
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Thus, if we can find a covering of the integers which uses the moduli 2 and 4 each 
once and the modulus 8 twice, we can obtain congruences on A as in Lemma 2.2 (ii) 
to guarantee that for a fixed choice of 6, some prime in {11, 73, 101, 137} divides each 
of the numbers A+ 6-10" for k € Z* U {0}. One checks that 


k 


1 (mod 2) 
k = 2 (mod 4) 
k = 0 (mod8) 
k = 4 (mod8) 


is such a covering of the integers. We take 6 = 7 and ensure Lemma 2.2 (ii) holds by 


taking 
A = -7-10' =7 (mod 11), 
A = -7-10? =7 (mod 101), 
A = -—7-10° =-7 (mod 73), 
A = -7-104 =7 (mod 137). 


Then since the conditions of Lemma 2.2 are satisfied, all of the numbers A+7-10* are 
divisible by some prime in the set P7 = {11, 73,101,137}. The Chinese Remainder 
Theorem implies that the four congruences above of the form A = —7- 10* together 
with the congruence A = —2 (mod 3) are all equivalent to a single congruence modulo 
the product of the primes in P2 UP; UP; U Pg = {3, 11,73, 101,137}. Taking A to 
satisfy this single congruence, we have that for 6 € {2,5,7,8} the numbers A+6-10*, 


where k € Zt U {0}, are divisible by some prime in 
P2U Ps UP, UP: = {3, 11, 73, 101, 137}. 


Furthermore, we have A is not divisible by primes in {3, 11, 73, 101, 137} so that there 
are primes p= A (mod 3- 11- 73-101 - 187). 
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We continue in this way, finding for each 6 € {1,3,4,6,9} coverings 
k=b; (mod c(p;)) 


with corresponding congruences A = —6 - 10% (mod p;). To allow these congruences 
to occur simultaneously with congruences already established, our choices will either 
utilize p; distinct from those used before or will be equivalent to congruences already 


established. For example, we have A = 7 (mod 11) above, so it follows that 
A=-—4-10° (mod 11). 


Recall c(11) = 2. Thus, taking k = 0 (mod 2) in (i) of Lemma 2.2 and taking 6 = 4, 
we see that (ii) of Lemma 2.2 holds for the prime 11. In other words, we may re-use 
c(11) = 2 as a modulus in a covering corresponding to 6 = 4 so long as we choose 
k =0 (mod 2). Table 2.1 exhibits such a covering, where in the table we have made 


use of the notation 
pio = 9999999900000001, pig = 15343168188889137818369, 


Ppoq = 915217525265213267447869906815873, —poq = 1253224535459902849, 
pos = 53763491189967221358575546107279034709697. 


To clarify, we are taking P, to be the set of 27 primes appearing in Table 2.1 
in the columns with heading “prime p;”. One checks that the columns with heading 
“congruence” in Table 2.1 describe a covering of the integers. Note that the least 
common multiple of the moduli appearing in this system of congruences is 768, so 
checking that a covering of the integers is indeed described in Table 2.1 amounts to 
verifying that each of the numbers 0,1,..., 767 satisfies some congruence in Table 2.1. 
One further checks that each prime number p; in the columns with heading “prime 
p;” indeed has c(p;) equal to the modulus appearing in that row. This suffices to 


show that the congruences in Table 2.1 satisfy (i) in Lemma 2.2. Thus, with 6 = 4 
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Table 2.1 Covering used in Lemma 2.2 (i) for A+ 4-10" 


F row congruence prime p; 
row congruence prime p; 
= 14 | 45 (mod 64) | 834427406578561 
i k =0 (mod 2) 11 
5 Eeiimnod 6) Z 15 | 15 (mod 96) 97 
aoe Sa Gadel) a 16 | 39 (mod 96) 206209 
7 i = 9 (mod 16) | 5882353 ii 63 (mod 96) 66554101249 
— 18 87 (mod 96) 75118313082913 
5 | k=3 (mod 32) | 353 
6 |k=11 (mod 32) | 449 
7 |k=19 (mod 32)| 641 ote AN ee 1) ue 
= 21 | 61 (mod 128) | 1265011073 
8 | k =27 (mod 32) 1409 
22 | 47 (mod 192) 193 
9 | k=5 (mod 32) | 69857 
= 23 | 95 (mod 192) 769 
10 | k = 23 (mod 48) P10 
Il |k=13 (mod 64) | 19841 
= 25 | 191 (mod 192) Dos 
12 | k =21 (mod 64) | 976193 
12 00 Gnoaoay WG1S7457 |. oo: | tee mod 206) aul 
— 27 | 253 (mod 256) 15361 


and taking A to be a solution to the corresponding congruences in Lemma 2.2 (ii), 


e.g., from row 9 we take 


A=-—4-10° (mod 69857), 


we have from Lemma 2.2 that for all k € Z* U {0}, the number A+4- 10* is divisible 
by some prime in Py. Thus, for any number A satisfying the system of congruences 
already established, we have for 6 € {2,4,5,7,8} that the numbers A+ 6-10" are 
divisible by some prime p € Pz U Py U Ps U Pz U Ps. 

For 6 € {1,6,9}, we utilize the remaining coverings appearing in [10], which we 
reproduce here in Tables 2.2, 2.3, and 2.4. The following notation for primes are used 


in these tables: 


p7 = 440334654777631, po = 5964848081, py, = 102598800232111471, 


Pig = 3199044596370769, pi3 = 265212793249617641, 


Pia = 30703738801, pis = 625437743071, 


Pig = 97802050308786191965409441, = py7 = 4185502830133110721, 


Le 


pa, = 4458192223320340849, 


p31 = 60368344121, 


Pa7 = 127522001020150503761, 


p34 = 73765755896403138401, 


p36 = 603812429055411913, 


p35 = 11189053009, 


p32 = 848654483879497562821, 


p37 = 148029423400750506553. 


Table 2.2 Covering used in Lemma 2.2 (i) for A + 1- 10* 


row congruence prime p; row congruence prime p; 
1 k =0 (mod 3) 37 8 | k =26 (mod 54) | 70541929 
2 k =1 (mod 6) 13 9 | k=53 (mod 54) | 14175966169 
3 k =2 (mod 9) 333667 10 | k=4 (mod 12) 9901 
4 | k =5 (mod 18) 19 11 | k = 10 (mod 24) | 99990001 
5 | k=14 (mod 18) | 52579 12 | k = 22 (mod 72) Pia 
6 | k=8 (mod 27) 757 13 | k = 46 (mod 72) 3169 
7 | k=17 (mod 27) P7 14 | k=70 (mod 72) 98641 
Table 2.3 Covering used in Lemma 2.2 (i) for A +6 - 10" 
row congruence prime p; 

5 Ee0 Gods) 7 row congruence prime p; 
5 i = T (mod 5) al 11 | k =18 (mod 60) 61 
Ee 20 (ed 1) 001 12 | k =38 (mod 60) | 4188901 
= 13 | k =58 (mod 60) | 39526741 
4 | k =3 (mod 20) 3541 — 
— 14 | k =4 (mod 15) 31 
5 | k=13 (mod 20) | 27961 = 
= 15 | k=9 (mod 15) | 2906161 
6 | k=7 (mod 30) 211 = 
= 16 | k=14 (mod 45) | 238681 
7 | k=17 (mod 30) 241 — 
8 | k=27 (mod 30) | 2161 phi Se no) 2s 
eae Gao 18 | k =44 (mod 90) | 29611 
= 9 = 
10 | E=28 (mod 40) Tavesr | 29 [F=89 (mod 90) | 2762001 


Taking A to simultaneously satisfy all the congruences corresponding to (ii) in 
Lemma 2.2, we have that for all 5 € {1, 2,4, 5, 6,7, 8,9}, each of the numbers A+6-10* 


is divisible by some prime p € P; UP2U P4U Ps UP. U Py U Ps U Po. We require one 


more covering system for 6 = 3. 


For 6 = 3, we will make use of three primes previously used for other choices 


of 6, namely the prime 73 for 6 = 7 and the primes 7 and 17 for 6 = 4. From 
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Table 2.4 Covering used in Lemma 2.2 (i) for A +9 - 10* 


row congruence prime p; ; 
; = 17 5a row congruence prime p; 
=0 (m 
su cunod. 20 | k=33 (mod 84) | 226549 
2 k =1 (mod 7) 4649 = 
— 21 | k=75 (mod 84) p21 
3 k =2 (mod 21) 43 = 
= 22 | k=19 (mod 63) | 10837 
4 k =9 (mod 21) 1933 
= 23 | k=40 (mod 63) | 23311 
5 k = 16 (mod 21) | 10838689 — 
= 24 | k=61 (mod 63) | 45613 
6 k =3 (mod 14) 909091 — 
= 25 | k=6 (mod 28) | 121499449 
7 k = 10 (mod 28) 29 = 
= 26 | k=13 (mod 56) 7841 
8 | k= 24 (mod 28) 281 = 
— at k= 41 (mod 56) P27 
9 k = 4 (mod 35) 71 — 
10 | k=11 (mod 35) | 123551 26) Idee and Eey) ae 
p= Ce 29 | k =48 (mod 140) | 3471301 
= ai 30 | k= 76 (mod 140) | 13489841 
12 | k=25 (mod 70) | 4147571 = 
13 | k=60 (mod 70) P13 — 
= 32 | k= 132 (mod 140) | pap 
14 | k =32 (mod 105) pia — 
= 33 | k =27 (mod 112) 113 
15 | k =67 (mod 105) P15 = 
16 | k =102 (mod 105) P16 — 
— 35 | k =55 (mod 168) P35 
17 k =5 (mod 42) 127 — 
18 | k&=26 (mod 42) | 2689 are = er Gad 16s) 
19 | k=12 (mod 42) | 459691 bls 


our prior congruences, we have A = —7 (mod 73), A = —4-10 = 2 (mod 7), and 


A =-—4-10=-—6 (mod 17). We make use of the equivalent congruences 


A = -3-10° (mod 73), 
A = -3-10* (mod 7), 
A = -3-10'° (mod 17) 


to help construct a covering for the case 6 = 3. Our choices for this covering will 
require primes where 10 has multiplicative order larger than that appearing in the 
tables in [3]. To find such primes, we obtain partial factorizations of ®,(10), where 
®,,(x) is the nth cyclotomic polynomial. We make use of the following facts regarding 
primes dividing ®,,(10), which are discussed in [3, p. III C 1]. Let p be a prime dividing 
®,,(10). 


e Ifp>n, we have c(p) =n and p=1 (mod n). 
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e Ifp<nandn > 2, then p is the largest prime dividing n and v,(®,(10)) = 1. 


To find a lower bound on the number of primes p having c(p) = n, we first determine 
whether the largest prime q dividing n also divides ®,,(10). We then look for distinct 
small primes p;,...,p, with p; = 1 (mod n) and p,|®,(10). Note that the condition 
p; = 1 (mod n) helps narrow the search space of primes to consider. Once this 
calculation is complete, we verify that gcd (®,,(10) V4 piss bi) = 1. We then use 
the is_prime_power() routine in Sage to verify that ®,,(10)/]]/_, p; (or alternatively 
®,,(10)/(qTTj_, pi) if g | ®,(10)) is neither 1 nor a prime power. If this quotient is 
1, we have at least r distinct primes p with c(p) = n. If the quotient is a prime 
power, we have at least r+ 1 distinct primes p with c(p) =n. If the quotient is not 
a prime power and not 1, we have at least r + 2 distinct primes p with c(p) = n. In 
the latter case, although the additional two primes can be determined with enough 
computation, we do not need them explicitly for the proof of Theorem 1.3. 

Table 2.5 displays the results of these computations of prime divisors of ®,,(10). 
Each n we consider has largest prime divisor g = 11 and only n = 242 has 11 | ®,,(10). 
In Table 2.5, we write ®©,(10) = pipo---p-C,, where C, is a composite factor of 
®,,(10) having at least 2 distinct prime divisors different from pj, po,...,p,. We 
write ®,(10) = pypo---prPn, where P, is a prime factor of ®,(10) different from 
D1, P2,-++,Pr- Computationally, P, was determined to be a prime power and then 
verified to be a prime. In the tables that follow, we denote by p,; and py» the two 
smallest prime divisors of C’,. We did not compute the values of pp, and pp», but 
we know they exist. We also included some values of n in Table 2.5 which appear in 
[3] since this notation allows us to avoid displaying lengthy primes. 

Table 2.6, Table 2.7, and Table 2.8 display the covering system of the integers 


corresponding to the case 6 = 3 in Lemma 2.2 (i). To conserve space, we set 


Pig = 16205834846012967584927082656402106953, 
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Table 2.5 Partial factorizations of ®,,(10) for large n 


n Partial factorization of ®,,(10) 
121 15973 - 38237 - 274187 - Pia) 
242 11 - 4357 - 25169 - 1485397 - Co49 
275 7151 - 15401 - 59951 - Coz5 
363 622001227 - C'363 
396 79082656489 - C396 
484 | 56629 - 170369 - 29606281 - 1491164086760128255001869 - C4g4 
605 C05 
726 727 - 1453 - C726 
792 761113 - Cr92 
1188 765144469 - Chigg 
1210 10891 - 131891 - Charo 
1452 Ps 
2420 1006721 - 2323201 - 1754328181 - C429 
2904 Cra 
4356 949609 - 384538969 - C4356 
5808 21582529 - 114690577 - Csgog 


pag = 138267770127916457629034873443951, 
pag = 170354891389249407509 7664562023844278044121, 
paa = 1395900370916327245555441901, 
pas = 36380545029953205956377406 702261, 
pag = 1112314101311286003379752617807870409611285281, 
P52 = 136614668576002329371496447555915740910181043, 
Pe = 362853724342990469324766235474268869786311886053883, 
pea = 141122524877886182282233539317796144938305111168717, 
p75 = 790700930759469400105355200058865839110097409345 7603716419437, 
P77 = 8927244623941181398233253, — p7g = 1866763546567680996103417376059, 
P79 = 112970308382439859401726947341740704554951737408511354573, 


ps0 = 7450755712209696453 10060665 14788984423438931950911932421947947339. 


pal 
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Table 2.6 First part of covering used for A +3- 10* 


row congruence prime p; row congruence prime p; 
1 k =4 (mod 6) 7 22 | k = 159 (mod 264) 6796152793 
2 | k=3 (mod 8) 73 23 | k = 247 (mod 264) 24387741577 
3 | k=0 (mod 11) 21649 24 | k=6 (mod 55) 1321 
4 | k=1 (mod 11) 513239 25 | k=17 (mod 55) 62921 
5 k = 10 (mod 16) 17 26 k = 28 (mod 55) 83251631 
6 k =2 (mod 22) 23 jas k = 39 (mod 55) | 1300635692678058358830121 
7 k = 13 (mod 22) 4093 28 | k =50 (mod 110) 331 
8 k =3 (mod 22) 8779 29 | k = 105 (mod 110) 5171 
9 k = 14 (mod 44) 89 30 k =7 (mod 110) 20163494891 
10 k = 36 (mod 44) 1052788969 31 | k=62 (mod 110) | 318727841165674579776721 
Il | k= 15 (mod 33) 67 32 | k= 18 (mod 220) 661 
12 k = 26 (mod 33) | 1344628210313298373 33 | k=73 (mod 220) 18041 
13 k = 37 (mod 66) 599144041 34 | k = 128 (mod 220) 148721 
14 | k = 38 (mod 132) 5419170769 35 | k = 183 (mod 220) 1121407321 
15 | k=126 (mod 132) 789390798020221 36 | k= 40 (mod 275) 7151 
16 | k =60 (mod 132) 2361000305507449 37 | k=95 (mod 275) 15401 
17 k =5 (mod 88) 617 38 | k = 150 (mod 275) 59951 
18 k = 49 (mod 88) Dig 39 | k = 205 (mod 275) P275,1 
19 | k = 104 (mod 264) 2377 40 | k = 260 (mod 275) P275,2 
20 | k = 236 (mod 264) 16369 41 | k =29 (mod 165) 471241 
21 | k=71 (mod 264) 432961 42 | k = 84 (mod 165) Par 
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Table 2.7 Second part of covering used for A +3 - 10° 


row congruence prime p; row congruence prime p; 
43 | k = 139 (mod 165) pag 64 | k=152 (mod 198) P64 

44 | k =51 (mod 220) Daa 65 k = 64 (mod 396) 79082656489 
45 k= 161 (mod 220) P45 66 k = 163 (mod 396) P396,1 

46 | k = 106 (mod 330) 4124507971 67 k = 262 (mod 396) P396,2 

47 | k = 216 (mod 330) | 19835636682880495867311241 68 k = 361 (mod 792) 761113 
48 | k = 326 (mod 330) Dag 69 k = 757 (mod 792) P792,1 

A9 k =8 (mod 77) 5237 70 k = 75 (mod 297) 55243 

50 k = 19 (mod 77) 42043 yl k =174 (mod 297) 198397 
51 k = 30 (mod 77) 29920507 ie k = 273 (mod 297) 1981560241 
52 k = 41 (mod 77) P52 73 k = 86 (mod 297) 31600574312077 
53 | k=52 (mod 154) 463 74 k = 185 (mod 297) | 165426670443186506567467 
54 | k = 129 (mod 154) 24179 75 | k =284 (mod 297) P75 

55 | k= 63 (mod 154) 590437 76 | k=97 (mod 594) 7129 

56 | k = 140 (mod 154) 7444361 77 | k=196 (mod 594) Piz 

57 | k=74 (mod 154) 4539402627853030477 78 k = 295 (mod 594) D78 

58 | k = 151 (mod 154) 4924630160315726207887 79 k = 394 (mod 594) D79 

59 k =9 (mod 99) 199 80 | k= 493 (mod 594) Pso 

60 k = 20 (mod 99) 397 81 | k =592 (mod 1188) 765144469 
61 | k=31 (mod 99) 34849 82 | k= 1186 (mod 1188) ete 

62 | k =42 (mod 99) Pe2 83 | k= 10 (mod 121) 15973 

63 | k=53 (mod 198) 7093127053 84 k = 21 (mod 121) 38237 
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Table 2.8 Third part of covering used for A +3 - 10* 


row congruence prime p; row congruence prime p; 
85 | k= 32 (mod 121) 274187 105 | & =593 (mod 1452) Pas 
86 | k = 43 (mod 121) Pai 106 | k = 1319 (mod 4356) | 949609 
87 | k=54 (mod 242) 4357 107 | k = 2771 (mod 4356) | 384538969 
88 [k= 175 (mod 242) 25169 108 | k = 4223 (mod 4356) | passoa 
89 | k = 65 (mod 242) 1485397 109 | k = 714 (mod 2904) | poooat 
90 | k = 186 (mod 242) P2421 110 | k = 1440 (mod 2904) P2904,2 
91 | k= 76 (mod 242) Doar,2 it | k = 2166 (mod 5808) | 21582529 
92 | k = 197 (mod 484) 56629 112 | k = 5070 (mod 5808) | 114690577 
93 | k =439 (mod 484) 170369 113 | & = 2892 (mod 5808) | pssoa 
94 | k = 87 (mod 484) 29606281 114 | k = 5796 (mod 5808) P5808,2 
95 | k = 208 (mod 484) | 1491164086760128255001869 115 | k=120 (mod 605) P605,1 
96 | & = 329 (mod 484) P4g4,1 116 | k = 241 (mod 605) P605,2 
97 | k = 450 (mod 484) Pag2 117 | k = 362 (mod 1210) 10891 
98 | k= 98 (mod 363) 622001227 118 | k = 967 (mod 1210) 131891 
99 k = 219 (mod 363) P363,1 119 k = 483 (mod 1210) P1210,1 
100 | k = 340 (mod 363) P363,2 120 | k = 1088 (mod 1210) P1210,2 
101 | k = 109 (mod 726) 727 121 | k = 604 (mod 2420) 1006721 
102 | k = 230 (mod 726) 1453 122 | k = 1209 (mod 2420) 2323201 
103 | k = 351 (mod 726) P726,1 123 | k = 1814 (mod 2420) | 1754328181 
104 | k = 472 (mod 726) P726,2 124 | k = 2419 (mod 2420) P2420,1 


Taking A to satisfy the corresponding congruences in Lemma 2.2 (ii) with 6 = 3, 
we have that A+3-10* is divisible by some prime p; appearing in the columns with 


) 


heading “prime p;” in Tables 2.6, 2.7, and 2.8. For the remainder of the paper, we 


will set 


M= [J IIp 


5€{1,2,...,8,9} pEPs 


and will take A to be the unique residue class modulo M satisfying all the congruences 
in Lemma 2.2 (ii) we have constructed for 6 € {1,2,...,8,9}. Note that each of the 
congruences A = —6 - 10% (mod p;) describe nonzero residue classes modulo p;. For 
p; = 3 and p; = 7, we avoided 6 = 0 (mod p;) and for all other p; appearing above, 
p; > 10 ensures that —d - 10% is nonzero. Thus, we have gcd(A, M) = 1. Based on 
our work so far, if p is a prime with p= A (mod M), then p +6 - 10* is divisible by 
a prime dividing M for every choice of 6 € {1,2,...,9} and k € Zt U {0}. We now 


proceed to 6 € {—9, —8,..., -2,—1}. 


2.2 DECREASING A DIGIT 


Let K be a large integer. Our goal is to estimate the number of digitally delicate 
primes p with p = A (mod M) and with p having < K digits. Observe that the 
number of digit changes for such p is at most 9K. As noted above, we have that 
pt+06-10', for0<j< K and6 € {1,2,...,8,9}, is divisible by a prime dividing M. 
For p > M and for 0 <j < K and 6 € {1,2,...,8,9}, we deduce then that p+6-10/ 
is composite. We extend this to the 9K more numbers p+ 6- 10% where 0 <j < K 
and 6 € {—9, —8,...,—-2, —L}. 

We begin with an idea of P. Erdés to dispense with most of these 9K possibilities 
by considering primes p only in residue classes modulo some primes for which 10 has 
small order. To be more precise, suppose kg is the smallest value of the order of 10 
modulo a prime divisor of M. Fix « > 0, and let Ky = Ko(e) be a positive integer 


to be specified momentarily. Note that with ko fixed as above, the value of Ko will 
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depend only on ¢, and K will be taken to be large in comparison to Ko. For each 
integer c € (ko, Ko], we know there exists a prime q. such that 10 has order c modulo 


de (cf. [1], [31]). A result of Stewart [29] implies that we can take 
de > ke: (log c)’, (2.1) 


for some constant k, > 0 (depending on ky). We proceed with the gq. so chosen. Since 
the order of 10 modulo primes dividing M are < ko, no q. divides M. 

We will restrict our attention to primes c € (ko, Ko]. We want to use most of 
these primes and separate them into 9 sets depending on 6 € {—9, —8,...,—2, —I}. 


For this purpose, we define 


Os = {c € (ko, Ko] : c prime,c = 6 (mod 11)}. 


We set 
Q= U Q5, M = II ike and W=MM"'. 
6€{-9,-8,...,-1} dE€{-9,-8,...,-1} cE Q5 
For each 6 € {—9, —8,..., —2, —1}, we consider 


p= -6 (mod |[f «). 


cE Qs 
By the Chinese Remainder Theorem, there is a B € Z such that these 9 congru- 
ences, one for each 6 € {—9,—8,...,—2,—1}, are equivalent to a single congru- 
ence p = B (mod M'). As a consequence, we have that if p = B (mod M’) and 
6 € {-9, -8,...,-2,—1}, then for every nonnegative integer k that is divisible by 
ac € Q;, we have p+ 6 -10* is divisible by the prime q.. Thus, with fixed p = B 
(mod M’) and fixed 6 € {—9, —8,..., —2, —1}, an inclusion-exclusion argument shows 


that the number of integers 7 € (0, K] with p+ 6-10’ relatively prime to M’ is 


<K I] ( = "| + 0(2%°), (2:9) 
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We take Ky = Ko(e) so that for K sufficiently large, the expression in (2.2) is << K, 
which is possible since }°?.,,, 1/c diverges where the * indicates the sum is over primes 
c= 0 (mod 11). 

Since gcd(M, M') = 1, the primes p which are A modulo M and B modulo M' 
correspond to the primes in a single congruence class, say w, modulo W = MM’. 
We now have that if p= w (mod W), then for all but < 9¢K choices of O< 7 < K 
and 6 € {—9,—-8,...,—-2,—-1} U {1,2,...,8,9}, the number p + 6 - 10” is divisible 
by a prime dividing W. Furthermore, the choices for 6 - 10’ for which p+ 6 - 10/ 
is divisible by a prime dividing W depend only on the residue class w modulo W 
and, hence, do not depend on the prime p = w (mod W). Let H denote the set of 
6-10? with 0 <j < K and 6 € {-9,-8,...,—-2,-1} U {1,2,...,8,9} such that the 
number p+ 6-10/ is not divisible by a prime dividing W for p= w (mod W). For 
our purposes, we can in fact take H to be the numbers of the form 6 - 10’, where 
6 € {-9, -8,...,-2,—-1} and j is not divisible by any c € Q5. Note that |H]|, the 
size of H, is bounded above by 9eKk. 

We continue with an idea of Tao [30] and recall a theorem from sieve theory (cf. 


[13, Theorem 2.4], [30]). 


Theorem 2.3. Let W and w be positive integers, and let h be a non-zero integer. 
If x is sufficiently large (depending on W and w), then the number of primes p < x 
with p= w (mod W) and p+h prime is 


Cz ie Vee 
< Fw )dogay? , LI ( 7 i) ? 


q\(hW) 


where C' is an absolute constant. 


Of some significance here is that C' does not depend on h. For our purposes, we want 
h = 6-10? € H. Thus, the primes dividing h belong to the finite set {2,3,5, 7}. 
We also view M and M’ and, hence, W as fixed. We take x = 10* with K large 


depending on ¢. For each of the h = 6-10? € H, we apply Theorem 2.3. 


yas 


Recall that the primes dividing W appearing in the product in Theorem 2.3 
include the primes q. for c € Q and Q depends on Ko and Ko depends on ¢. We 
will use a lower bound on this product appearing in the theorem to get rid of its 
dependence on ¢. More precisely, recalling (2.1) and that the primes dividing h 


belong to the set {2,3,5,7}, we see that 


1 8 1 8 1 1 
Pee | Se i een nee pe Nees 
Tr ( ;) ~ 39 TL ) 35 an ;) H ( *) 


q|(hW) q\W p|M peek Ky 
8 1 1 
>= 1-—- eS eS 
395 ae ( *) role ( kyc(log =) 
p\|\M 


As the sum ¥>,5,;, 1/(c(logc)”) converges, we see that, for a fixed h € H, there 
is a constant C’ = C’(M), independent of « > 0, such that there are at most 
C'x/(¢(W) log? x) primes p < x with p = w (mod W) and p+h prime. Letting 
h vary over the elements of H, the total number of primes p < x with p = w 
(mod W) and p+h prime for some h € H is bounded above by 


|H|- C's z eK Cre: 9eC"x 
o(W) log? x ~ ¢(W) log? —(W) log(10) - log x’ 


The remaining primes p < x with p = w (mod W) satisfy that p + h is composite 
for every h € H. Thus, these remaining primes will be digitally delicate. By the 
Prime Number Theorem for Arithmetic Progressions, the total number of primes 
p <« =10* with p = w (mod W) and K large is > x/(2¢(W) log). By taking 
€ sufficiently small and, hence, K sufficiently large, we see that there are a positive 
proportion of primes p < x = 10" with p= w (mod W) that are digitally delicate. 
As these primes satisfy p = A (mod M), we deduce that a positive proportion of 
primes up to 10" are widely digitally delicate. As a positive proportion of primes 
up to 10* is still a positive proportion of primes up to 10**+!, we can see that the 
restriction of x to the form 10“ can be eliminated and there is a uniform lower bound 
on the proportion of primes up to x, for x large, that are widely digitally delicate. 


This establishes Theorem 1.3. 
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2.3. RELATED TOPICS AND OPEN PROBLEMS 


Questions related to Theorem 1.3 abound. Carl Pomerance (private communication) 
has asked about primes which are not digitally or widely digitally delicate. In this 
regard the simplest question one could ask is whether there are infinitely many primes 
which are not digitally delicate. The prime k-tuple conjecture states that every 
admissible pattern for a prime constellation occurs infinitely often. If one assumes 
this conjecture, then in the case k = 3 we obtain there are an infinite number of 
primes p for which p — 2 and p+ 4 are also prime. Then for any base b > 6 express 
such p in base 6 as 


where 0 < d; < b—1. Then if do > 2, we have that p — 2 corresponds to changing 
a digit of p in base 6. If alternatively dg < b — 4, we obtain that p + 4 corresponds 
to changing a digit of p in base 6. Since for b > 6 we must have either dp > 2 or 
do < b—4, the prime p in such a constellation cannot be digitally delicate in base b. 
Thus an infinite number of primes are not digitally delicate. 

Apart from utilizing a conjecture, we have been unable to establish an uncondi- 
tional proof of the infinitude of primes which are not digitally delicate. What is sur- 
prising in this regard is that computationally it appears as if the proportion of digitally 
delicate primes is less than the proportion of non-digitally delicate primes. For each 
of the bases b € {2,3,--- ,10} and each of the intervals (10!°°, 101°), (102°, 101°°), 
(102°, 1017), (1017, 107°), we generated 100000 primes in the given interval and 
counted the number of primes out of the 100000 that were digitally delicate in that 
base. We use the notation L,(k, 2) in Table 2.9 to represent the number of primes 
out of 100000 in the interval (10*, 10°) that were digitally delicate in base b. 

The proportions seem to be constant for a fixed base, with the largest proportion 


appearing when b = 3 at around 21% of primes being digitally delicate in base 3. 
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Table 2.9 L,(k, 2) for bases b € {2,3,--- ,10} 


b 2 3 4 5 6 c 8 9 10 
L£,(100, 125) | 16246 | 20860 | 4159 | 18166 | 75 | 9976 | 375 | 1520 | 18 
£y(125, 150) | 16319 | 20641 | 4285 | 18171 | 68 | 10270 | 406 | 1513 | 12 
£y(150,175) | 16352 | 20736 | 4383 | 18396 | 66 | 10063 | 417 | 1510 | 15 
£y(175, 200) | 16385 | 20747 | 4280 | 18212 | 65 | 10223 | 376 | 1474 | 12 


However, perhaps this pattern changes as the primes involved grow. Such a determi- 
nation could be made with more computational effort. 

One discernible pattern in Table 2.9 that bears out theoretically is the diminishing 
proportion of primes that are digitally delicate in bases which are non-prime powers 
of primes. One can see in Table 2.9 that the proportions for 4, 8, and 9 are much 
smaller than the proportions for any of the primes 2, 3, 5, 7. This is because a change 
of a single digit of p in base b implies a change of a single digit of p in base b’, so that 


the following holds. 


Lemma 2.4. Let b be an integer greater than 1 and let p be a prime digitally delicate 


in base b', t>1. Then p is digitally delicate in base b as well. 


This result does not hold for arbitrary composite bases. For example, the prime 
28151 is digitally delicate in base 6 but not in base 2 or base 3. Despite this non- 
implication, it appears from the columns b = 6, b = 10 in Table 2.9 that the proportion 
of digitally delicate primes decreases as the number of divisors of b increases. 

Moving beyond the basic question of their infinitude, one could also ask whether 
a positive proportion of primes are not digitally or not widely digitally delicate. Since 
the set of primes not digitally delicate is contained within the set of primes not widely 
digitally delicate, perhaps it is easier to show that there are an infinite number of 
primes not widely digitally delicate, or even a positive proportion of such primes. 
However, we have not made progress in this regard, and expect that methods other 


than those used in this dissertation would be necessary to establish such a result. 
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In [10], the existence of infinitely many digitally durable composite numbers n 
was established, that is, composite n satisfying the property that if any digit of n 
is replaced by a different digit, then the resulting number is composite. This easily 
follows from the methods of Erdés [5]. Extending this idea to the leading 0’s of n, 


we define widely digitally durable composite numbers in the natural way. 


Definition 2.5. A composite number n is widely digitally durable if n has the prop- 
erty that if any digit of n, including any of its infinitely many leading 0’s, is replaced 


by a different digit, then the resulting number is composite. 


It is not difficult to show that the approach in this dissertation also leads to the 


following result regarding such composite numbers. 


Theorem 2.6. There is a positive proportion of composite numbers that are widely 


digitally durable in base 10. 


Proof. Take A (mod M) as in the proof of Theorem 1.3 and pick any prime r not 


dividing M. Let A’ be the unique solution modulo r - M to the system 
x=A (mod M), x=0 (modr). 


Thus any number m satisfying m = A’ (mod r- M) has both the property that m is 
composite as well as the property that increasing any digit of m, including any of its 
infinitely many leading 0’s, results in a number divisible by some prime dividing M. 

The proof proceeds similarly to that in Section 2.2. Let K be a large integer. 
Our goal is to estimate the number of digitally durable composites d with d = A’ 
(mod rM) and with d having < K digits. Observe that the number of digit changes 
for such d is at most 9K. We have that d+ 6-10’, for 0 < j < K andoé € 
{0,1,2,...,8,9}, is divisible by a prime dividing rM. Ford >rM and for0 <j < K 
and 6 € {0,1,2,...,8,9}, we deduce then that d+6-10/ is composite. We extend this 


to the 9K more numbers d+ 06-10% where 0 < j < K and 6 € {—9, —-8,...,-2,—1}. 
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We continue as in Section 2.2 with the ideas of P. Erdos to dispense with most 
of these 9K possibilities by considering composites d only in residue classes modulo 
some primes for which 10 has small order. To be more precise, suppose ko is the 
smallest value of the order of 10 modulo a prime divisor of rM. Fix ¢ > 0, and let 
Ko = Ko(e) be a positive integer to be specified later. For each integer c € (ko, Kol, 
we choose the prime q, as in Section 2.2 so that 10 has order c modulo q,. Since the 
order of 10 modulo primes dividing rM are < ko, no q. divides rM. We restrict our 
attention to primes c € (ko, Ko] and again define for 6 € {—9,—8,--- ,—2,—1} the 
set 


O5 = {c € (ko, Ko] : c prime,c = 6 (mod 11)}. 


Fix, as before, 


Q= U 95, MS II II « and W=MM"’. 
6€{-9,-8,...,-1} bE€{-9,-8,...,-1} cE Qs 
For each 6 € {—9, —8,..., —2, —1}, we consider 


d=-—6 (mod [I «). 


c€O5 


By the Chinese Remainder Theorem, there is a B € Z such that these 9 congru- 
ences, one for each 6 € {—9,—8,...,—2,—1}, are equivalent to a single congru- 
ence d = B (mod M'). As a consequence, we have that if d = B (mod M') and 
6 € {-9, -8,...,-2,—1}, then for every nonnegative integer k that is divisible by 
ac € Q;, we have d+ 6 -10* is divisible by the prime q.. Thus, with fixed d = B 
(mod M’) and fixed 6 € {—9, —8,..., —2, —1}, an inclusion-exclusion argument shows 
that the number of integers 7 € (0, K] with d+6-10 relatively prime to M’ is 


<K ITI ( = | +0(2*°), (2.3) 


cE Os 


We take Ko = Ko(e) so that for K sufficiently large, the expression in (2.3) is << «Kk, 
which is possible since }°%.,,, 1/c diverges where the * indicates the sum is over primes 


c=0 (mod 11). 
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Since gcd(rM, M’) = 1, the numbers d which are A’ modulo rM and B modulo 
M' correspond to the primes in a single congruence class, say w, modulo W =rMM’. 
We now have that if d= w (mod W), then for all but < 9eK choices of O< 7 < K 
and 6 € {—9, -8,..., -2, -1}U{0,1,2,...,8,9}, the number d+6-10/ is divisible by a 
prime dividing W. Furthermore, the < 9¢K choices for 6-10’ for which d+6-10/ is not 
divisible by a prime dividing W are such that 6 € {—9, —8,...,—2, —1} and j is free of 
prime factors from Q;. Let H denote this set of 6-10’ with 6 € {—9, —8,...,-2, —1} 
and integers 7 € (0, A’) free of prime factors from Q5. Note that |H]|, the size of H, 
is bounded above by 9eK. 

Now there are at least 10% /W — (10%) composite numbers d up to 10* satisfying 
d=w (mod W). For each of the numbers h € H, the number of composite d so that 
d +h is prime is at most 7(10"). Thus the number of composites d = w (mod W) 


so that d+6-10* is composite for all 5-10" = h € H is at least 
10* 10* 10* 1.2762 
—— — (9K + 1)r(10***) > —— — (9K +1 
qo Oe Sp SE) co) ( aa 

10% 2(9ek +1)10*% — 10% NG 

go es aie 

W K log(10) Ww K 
Taking € < 1/(18W), we then obtain that the number of such composites d is at least 
1Oe OP SO 
Ww 2W KO 


> 


Taking K sufficiently large, we see that there are a positive proportion of composites 
d <a =10* with d= w (mod W) that are digitally durable. As these composites 
satisfy d = A (mod M), we deduce that a positive proportion of composites up to 
10* are widely digitally durable. As a positive proportion of composites up to 10* is 
still a positive proportion of composites up to 10+, we can see that the restriction 


of x to the form 10" can be eliminated and there is a uniform lower bound on the 


proportion of composites up to x, for x large, that are widely digitally durable. 


The proof of the analogue to Theorem 2.6 follows similarly in bases b where 


suitable coverings and congruence classes A (mod M) have been found that allow for 
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increasing the digits of a number, but these methods do not hold for an arbitrary 
base. 

Beyond the question of digit changing, one can also ask instead about inserting a 
digit. In [10], infinitely many composite numbers n were shown to have the property 
that if an arbitrary digit is inserted in n or before or after n, then the resulting number 
is also composite. This leads to several questions. Are there infinitely many primes 
p or a positive proportion of primes p such that if an arbitrary digit is inserted in p 
or before or after p, then the resulting number is composite? A list of such primes 
in base 10 is accessible through [28], but it is not known if the list is infinite. Do a 
positive proportion of primes have this property? What about the analogous question 
with leading 0’s? We do not know the answers to these questions and expect that 
methods different than those used in this dissertation are required to establish such 


a result. 
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CHAPTER 3 
IRREDUCIBILITY CRITERIA BASED ON DEGREE FOR 


POLYNOMIALS WITH NON-NEGATIVE COEFFICIENTS 


3.1 PRELIMINARY RESULTS 


We begin with a few results that both motivate the steps we take and help establish 


a useful result. 


Lemma 3.1. Fiz an integerb > 1. Let f(x) be a polynomial with non-negative integer 
coefficients such that f(b) is prime. If f(a) is reducible, then f(x) has a non-real root 
in the disc Dy = {z € C: |b—2| < 1}. 


Proof. Assume that f(a) is reducible. Then f(x) = g(x)h(x) where g(x) and h(x) 
are in Z|x], have positive leading coefficients, and are not identically +1. Since f(b) 
is prime, we may take, without loss of generality, g(b) = +1 and h(b) = +p. Let c be 
the leading coefficient of g(x), and denote its roots including multiplicity (1, ..., 6,. 


Thus, the degree of g(x) is r and we have 


1 = |9(b)| = lel [] Jo - 85| = [] | - Gil. 
jel j=l 


Thus, we conclude that at least one of the roots of g(x), and hence of f(x), is in Dy. 


The lemma then follows when we recall that f(2) has no positive real roots since it 


has non-negative coefficients. 


One example of the usefulness of Lemma 3.1 can be found in the following lemma, 


which we use extensively in the proof of Theorem 1.4. 
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Lemma 3.2. Let n be a positive integer. A complex number a = re” with r > 0 
and0 <6 <12/n cannot be a root of a non-zero polynomial with non-negative integer 


coefficients and degree < n. 


Proof. Let f(x) = Xo a,x) be a non-zero polynomial in Z[x] of degree s with each 
a; > 0 and s <n. By way of contradiction, assume a = re” is a root of f(x) with 
r >0Oand 0 < 6 < a/n. Observe that a has a positive imaginary part. Since the 
complex conjugate of a is also a root of f(x), we see that s > 2. For 1 <k <n, we 


have 0 < k@ < 7, so 
Im(a*) =r*sin(k0)>0 for 1<k<n. 


Since f(x) has non-negative coefficients with deg f = s < n, and a has positive 
imaginary part, we have 


Im(f(a)) 2 Im(a*) > 0, 


contradicting the fact that a is a root of f(z). 


Having established Lemmas 3.1 and 3.2, a motivating idea throughout the rest 
of this chapter is to replace the disk 9, in Lemma 3.1 with a different region such 
that if a = re” is in this region, then |@| is bounded above by a small number. 
This combined with Lemma 3.2 will give sharp bounds, D,(b), such that if f(a) isa 
polynomial with non-negative integer coefficients with f(b) prime for integers b > 2 
and if the degree of f(x) is < D,(b), then f(x) must be irreducible, as discussed in 


Section 3.3. 


3.2 A ROOT BOUNDING FUNCTION 


For a given integer b > 5 and n € {3,4,6}, our main goal is to establish the bounds 
D,,(b) in Theorem 1.4. In this section we introduce certain rational functions that 


will give us information about the location of possible roots of f(x) assuming f() is 
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reducible. While better rational functions can be chosen for small b, as in [4], we will 
make choices to simplify the results for the purposes of this dissertation. 
Recall that ®,,(2) denotes the nth cyclotomic polynomial, and let ¢, = e?™/”. 
Fix an integer b > 5, and let f(x) be a non-constant polynomial with non-negative 
integer coefficients such that f(b) is prime. Suppose f(a) = g(x)h(x), where g(x) and 
h(x) are in Z[x], have positive leading coefficients, and are not identically +1. Since 
f(b) is prime, we may take, without loss of generality, g(b) = +1 and h(b) = +f(0). 
Using the ideas of [9], we want to show that either g(x) has a root in common with 


one of 
@3(¢ —b) = 2? — (2b-1)e +B? —b4+1, 


@,(2 — 6) = 2? — 2b + +1, 


(7 — b) =a" — (264+ le +6? +641, 


or g(x) has roots in a certain region R, to be defined shortly. 


We define 


F,(z) = (3.1) 
where 
Ns(z) = (lb + G3 — 2] + G3 — 2))?(lb +4 — 2||b — 4 — 2|)*([b + Ge — 2l|b + Ge — 21)? 


and 


D,(z) = |b—2z|*®. 


In the notation of [4], this amounts to choosing 
(€2(b), €3(0), ea(), e6(b), d(0)) = (0, 1,1, 1,1), 


but since our choice is constant we have simplified the presentation somewhat. 
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Setting z = x + iy, direct computations show that the expressions in Mj, and Dy 


simplify to 


(Jb + ¢3 — 2||b + G — 2|)? = y+ (2(@ — b)? + 2(@ — b) — 1)y” 


+ ((«— 6)? +2-—6+1)?, 
(Jb +4 —2||b—a— 2|)? = y+ (2(@ — 6)? — 2)y’ + ((@ — 8)? +1), 


(Jb + G5 — z||b + Go — 2|)” = y+ (2(a — 6)? — 2(@ — b) — Ly? 


+ ((« —b)? — (w— 6) +1)’, 
and 
lb—2/? a+ (a — 6). 


Notice that each one of these expressions is in Z[b, x, y?]. Thus, \,(z) and D;(z) are 
in Z[b, x, y*], making F;,(z) a rational function in b, 2, and y*. Moreover, we observe 


that for each integer b > 2, the polynomial 


Pr(x,y) = Dela + iy) — No(x + ty) (3.2) 
can be written as 
8 « 
Pi(x,y) = S> aj(b, ay” (3.3) 
j=0 


where each a;(b, x) is an integer polynomial in b and «. We write g(x) in the form 
g(x) = |] (x — 4), 
j=l 
where c is the leading coefficient of g(x) and (1,...,,, are the roots of g(x). Since 
we have f(x) = g(x)h(x), the 6; are also roots of f(x). It can be shown, as is shown 


in more generality in [4], that 


|g(b + ¢3)9(b + C3) 71g (b + a) g(b — t)|?|9(b + Go) g(b + Co)? 
Ig(b)|*° 
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and 
5 IL Ft) 
ct 
are equal. We denote this common value by V = V,(g). 

Since each of g(b+¢3)g(b+¢3), g(b+i)g(b—1), and g(b+¢6)g(b+,) is a symmetric 
polynomial in the roots of an irreducible monic quadratic in Z|x], we conclude that 
each of these expressions are themselves integers. Furthermore, g(b) = +1. Thus, 
by looking at the first expression for V, either V = 0 or V € Z*. In the case that 
V =0, we see from either expression for V that g(b + ¢3)9(b + C3), g(b + i)g(b — 1), 
or g(b + )g(b+ Gg) is zero. This happens precisely when g(x) is divisible by at least 
one of &3(x% — b), ®4(a — b), or B(x — b), so that if one of these shifted cyclotomic 
polynomials is not a factor of g(a), we have V € Zt. Observe that F;(z) is a non- 
negative real number for all z € C. By looking at the product in the second expression 
for V, we see that if V 4 0, then F,((;) > 1 for at least one value of j € {1,...,m}. 
Said in another way, if V #0, then there is a root §; of g(a), and thus of f(x), such 
that F,(8;) > 1. 

Summarizing the above ideas, given only that g(x) € Z[z], g(x) # +1, and 
g(b) = +1, we have shown that either g(x) is divisible by at least one of ®3(x — b), 


®,(x — b), and ®¢(x — b), or g(x) has a root £ in the region 


In the latter case, we will use an analysis of the region R, in the complex plane 
to obtain important information about the location of 6. For an illustration of Rp, 
Figure 3.1 depicts the region Rs. 

While analyzing the region R,, we will sometimes refer to points (x,y) in Ry. 
This is to be interpreted as the point z = x + zy in the complex plane in Ry. To 


further help analyze Ry», we consider P,(x,y) defined in equations (3.2) and (3.3). 
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The definition of D,(z) implies that D,(z) > 0 for all complex z 4 b. Thus 
Fi(a+iy)>1 and P,(z,y) <0 


are equivalent for z 4 b, and F,(x+iy) = 1 and P;(x, y) = 0 are equivalent for z # b. 
Note that P,(b,0) = Dy(b) — N4(b) = 0 — 1 = —1. Therefore, the z € C such that 
F(z) = 1 correspond exactly to the points (x,y) where P,(z, y) = 0. 


The following lemma corresponds to {9, Lemma 2], and the variation [4, Lemma 


Sek 
Lemma 3.3. Fix an integer b > 2. Then there exist real numbers ag = ao(b), 
a, = a,(b), and a non-negative real-valued function pp(x) defined on the interval 


Ty = |b — ao,b + a4] such that: 


oe 


(i) Po(x,y) £0 for alla ¢€ I, andy € 


(ii) Py(x, po(x)) = 0 for alla € Ih. 
(iii) pp(b — ap) = 0 and m(b+ a1) =0. 


(iv) The function pp(x) is continuously differentiable on the interior of I, and is 


continuous on Ip. 
(v) [fx andy are real numbers for which P(x, y) <0, then x € I, and |y| < py(x). 


In view of the above lemma, complex numbers of the form x+i;(x) are boundary 
points of ,, which are on or above the real axis. Since P;(x, y) is a polynomial in 
y” with coefficients in Z[b,z], the region R, is symmetric about the real axis. Thus, 
the points x — ip,(x) are boundary points of , which lie on or below the real axis. 
The points 6 — ao and b+ a; are boundary points on the real axis. 

A proof of Lemma 3.3 in more generality can be found in [4, Lemma 3.1]. In that 
proof, Cole et al. reduce the dependence on 6 to strictly a dependence on previously 


defined exponents €2(b), e3(b), e4(b), e6(b), and d(b). They then use Sturm sequences 
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on the resulting polynomials in x to prove the lemma. One can verify computationally 
that the desired properties hold for our choice of (€2, €3, e4, e6, 2) = (0,1, 1,1,1) when 
6 > 5, with ap = 1.522009... and a, = 1.522009.... Thus, the lemma follows. 

In the next section we will use Lemma 3.3 to prove irreducibility criteria based 


on the degree of f(z). 


3.3. IRREDUCIBILITY CRITERIA BASED ON DEGREE 
For the remainder of the chapter, we set for n € {3,4,6} and b € Z* the values 


6, =On(b) =arg(b+C,) and Dy = Da(b) = =| 


Because D,, features prominently in what follows, we explain here how the following 


lemma is a consequence of [20, Corollary 3.12] which states that the only rational 


values of cos(ar) with r € Q are 0, +1/2, and +1. 


Lemma 3.4. Let b € Z be greater than 2 and let n € {3,4,6}. Then 


T 
— ¢Z. 
at 


Proof. Letting r = arg(b+ ¢,)/a, it suffices to show r ¢ Q. Observe that for n € 
{3,4,6} we have 


arg(b+¢,) =arctan(z) where xe 


V3 1 v3 
2b-—1'b 2b+ 1] 


In each case, we have sin?(arctan(x)) = x?/(x* + 1) € Q, and hence 
cos(27r) = cos(2arg(b+ ¢,)) =1—2sin?(arg(b+¢,)) € Q. 


By Corollary 3.12 in [20], we deduce that if r € Q then cos(2 arg(b+¢,)) is an element 
of {0,+1/2,+1}. One checks that arg(3 + ¢,) < 7/6 for each n € {3,4,6}. Since 


arg(b+¢,,) decreases to 0 as b increases, we see that cos(2 arg(b+¢,)) ¢ {0, 41/2, +1} 


and the lemma follows. 


Al 


From Lemma 3.4, we obtain 


8 1 
i, | is arg(b+ Cn) 


One goal of this section involving D,, is to establish the following theorem. 


for b > 2 and n € {3,4, 6}. (3.5) 


r=| 


Theorem 3.5. Let b > 5 be an integer and Dy = |n/arg(b+G,)|. Let f(x) 
be a polynomial with non-negative integer coefficients such that f(b) is prime. If 


deg f(a) < Du, then f(x) is irreducible. 


While we will establish several generalizations of Theorem 3.5, the principle gen- 
eralization being Theorem 1.4, we begin with Theorem 3.5 because it is the product 
of a lengthy history in the literature. Its original predecessor appears as part of |7, 
Theorem 8], where Filaseta showed that the result above holds for all b > 2 when Dy is 
replaced by the value |7/arcsin(1/b)|. The values D4(b) and |7/arcsin(1/b) | usually 
coincide, but for some 8, e.g., values b that appear as the denominator of an odd con- 
vergent to the simple continued fraction for 7, we have D4(b) = |a/arcsin(1/b)| +1. 
Thus Theorem 3.5 is a sharpening of that aspect of [7, Theorem 8]. Later, in [4], 
Theorem 3.5 was established for integers 2 < b < 20, so we need only consider b > 20. 


However, the methods discussed here apply to all b > 5. 


Proof. We begin by showing that if deg f(x) < D4 in Theorem 3.5, then f(x) is 
irreducible. Recall the setup we have been using. That is, suppose that f(x) is 
reducible and write f(x) = g(x)h(x) where both g(x) and A(z) are in Zz], g(x) £1, 
h(a) # +1, and both g(x) and h(x) have positive leading coefficients. Furthermore, 
without loss of generality, we suppose that g(b) = +1. In Section 3, we established 
that either g(x) has a root in common with at least one of @3(x — b), ®4(a — 6), and 
¢(x — b), or g(x) has a root B = re” € Ry, where R, is defined as in (3.4). We also 
define P, (x,y) as in (3.2). 

We continue in two steps. First, we show that g(x) cannot share a root in common 


with any one of ®3(a — b), ®4(a — b), and ®¢(x — b) using Lemma 3.2. Second, we 
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will bound |6| using Lemma 3.3, then use Lemma 3.2 to show that g(x) cannot have 
a root in Rp, giving a contradiction. 

For the first step, suppose g(x) has a root 6 = re” = b+, where n € {3, 4, 6}. 
Then we have 0 = 6,, = arg(b+ ¢,,). To apply Lemma 3.2, we prove the following 


lemma. 


Lemma 3.6. Let b > 4 be an integer and D4, be as in Theorem 3.5. Then 
arg (b+ Gn) < a/D4 
for alln EN. 
Proof. For n € {1,2} and n > 4, the inequality 
arg (b+ Gn) < arg (b+ Gy) (3.6) 


is easily verified for all b > 4. For n = 3, we show that (3.6) also holds for b > 4. 


Observe that 


vV3/2 
ayy) 


arg (b+ ¢3) = arctan ( 5 


and arg (b+ ¢4) = arctan (=) 


Thus, for n = 3 and b > 4, we want 


V3/2 1 
<: ae 
arctan E = 1/2 < arctan (5) : 


which holds if and only if 


The latter inequality holds for b > 3.74. Thus, (3.6) holds for all b > 4 and n > 1. 


So for b > 4 and all n € N, we deduce that 


T T 7 


Fer Caer i - oe 


arg (b+ Cn) < arg (b+ Ca) = [m/arg(b+¢4)] Da’ 


where the strict inequality follows from (3.5). Thus the lemma follows. 
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By Lemma 3.2, since 0 < arg(b+¢,) < a/D4, any polynomial with non-negative 
integer coefficients with degree < Dy, cannot have a root whose angle is arg(b + ¢,). 
Thus, with f(x) as in Theorem 3.5 of degree < Dy, neither f(x) nor the factor g(x) of 
f(x), can have a root in common with any one of ®3(x — b), ®4(a —b), and ®¢(x—b). 
This completes the first step. 

For the second step, we consider the region R, and suppose that g(x) has a root 
B=re" € Ry, where we suppose as we may that r > 0 and 6 € (0,7/2). Let Ly be 
the line that passes through b+ ¢4 and the origin given by y = z/b. Figure 3.1 shows 
Rs with the line L;. Observe that, by definition, for all b > 5, the region FR, will be 
the same region as R; but shifted to be centered at (b,0). We claim that the line Ly, 
defined by y = x/b, lies completely above the region R, for b > 5. We use a Sturm 
sequence to see that Ps(z,2/5) has no real roots. Using Lemma 3.3 (ii) and (iv), we 
deduce that the line Ls; does not intersect the region Rs. Using a Sturm sequence 
with Ps(x,1), we deduce similarly that the line y = 1 does not intersect the region 
Rs. Since Ry is simply a horizontal translation of Rs, we see that the line y = 1 lies 
strictly above R, for b > 5. The lines y = 1 and Jy intersect at b+ C4, that is, the 
point (b,1). The smaller angle between these two lines is arctan(1/b) < arctan(1/5). 
It follows that, since the line Ls lies strictly above the region Rs, we have also that 
the line Ly lies strictly above the region Rp. 


Since 6 = re” € Ry with @ € (0, 7/2), we have that 
0<@<arg(b+G) <a/Dy (3.8) 


where the last inequality holds from Lemma 3.6. Thus by Lemma 3.2, f(x), and 


therefore g(x), cannot have a root in the region R,. This completes the second step 


and thus the proof of Theorem 3.5. 


We defer a discussion of the sharpness of the value D, until later, after we have 


established Theorem 1.4. The bound D, in Theorem 3.5 can be improved via The- 
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Figure 3.1 Ls; together with Rs 


orem 1.4 if one modifies the conclusion to allow for divisibility by one of the three 
cyclotomic polynomials @3(a — b), ®4(x — b), and ®¢(a — b). To prove Theorem 1.4, 
we first show the following lemma. 
Lemma 3.7. Let b be an integer > 5. Then 
7 us T 
D4g(b) < —~ < Dz (b) < —~ < De(b) < —> 
4 ( ) 6,(b) at ) 03(b) 6 ( ) 06(b) 


Proof. Note that for n € {3,4,6} and b > 5, the inequalities 


for b>5. (3.9) 


follow from (3.5). To establish the other inequalities, we write 6, for 0,,(b) and define 


mw _-mRe(b+ Ga) 7 
£n aera Im (b+ ¢,,) . 
and 
4, (b) = T 1 


~ 6, arctan (Im (b+¢,) /Re(b+ G,))’ 


Recall that 
Re(b+¢,) =b+cos(2a/n) and Im(b+¢,) =sin(27/n). 
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Thus, one can verify that £, (b) and LU, (b) are real, differentiable functions of b for 


b> 5. Using (3.5), and since x < tan(x) for x € (0,7/2), we have 
£y (0) < Dp (6) < Sty (6) for b€ Zand b> 5. (3.10) 
Next for n,m € {3, 4,6} define 
Brym (b) = Ln (b) — Un (8). 


Then we obtain 


OS ire 7 m Im (b a Ga 


db Im (b-+.G,) 7 lb+ Gale arctan (Im (b+ Gn) / Re (b+ Gaye 


Here, to simplify this derivative, and later, we use the Shafer-Fink ({11], [25]) inequal- 
ities 


3 
S a for x > 0. (3.11) 


—______ «< arct ae 
1+2V1+4+ 2? arcane) 14+2V1+4+ 2? 
For this proof, we will only need the lower bound due to Shafer [25]. We deduce that 
2 
dBm r rim (b+6m) (1+ 2y1+ (im (6+ Gm) /Re (b+ Gm))” 
db Im (6 + ¢,) b+ 6m 3Im (b+ Gn) / Re (b+ Gm) 
1 Tm+Gn) (Rel Cm) +2 |b , 


3Im (b + Gm) 


T mm tee 3 |b + Gn| ) 
Im (6 + Cn) Bae 31m (b + Gm) 


~Im(b+G,):  Im(Oo+c,) 


1 1 
cal (= (Q27/n) sin aa) 
One then verifies from the above that ¥3,4(b) and §¢3 (b) are increasing functions 


of b for b > 6. Since 3,4 (7) > 0.4 and 63 (7) > 2.4, for integers b > 6 we have 


One checks that (3.9) holds when b = 5 and b = 6. Thus, combining (3.5), (3.10), 


and (3.12), we deduce that (3.9) holds. 
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We now improve Theorem 3.5 with Theorem 1.4 by allowing for f(z) to be divisible 
by one or more of the cyclotomic polynomials ®3(x — b), ®4(a — b), or ®g(x — b). We 


restate Theorem 1.4 here for reference. 


Theorem 1.4. Fix an integer b > 5, and for n € {3,4,6} set 


T 


B= Dib" laepeey] m4 P= P= an =) | 


1000(25+1) 


Let f(x) € Z|a] with non-negative coefficients and with f(b) prime. If the degree of 
f(x) is < Dg, then f(x) is irreducible. Additionally, if the degree of f(x) is < D3 
and f(x) is reducible, then f(x) is divisible by ®4(a—b) and not divisible by ®3(x—b) 
or g(x — b). Furthermore, in the case that b > 26, if the degree of f(x) is < Dg and 
f(x) is reducible, then f(x) is divisible by either ®4(x — b) or ®3(x — b) and not by 
Oo(x — b). Lastly, in the case that b > 26, if the degree of f(x) is < D and f(x) is 
reducible, then f(x) is divisible by ®4(x — b), ®3(x — b), or ®3(ax — b). 


Before proceeding to the proof, we note that Cole et al. have already established 
the portions of Theorem 1.4 that hold for 5 < b < 20 as well as an analogous statement 
for D4 regarding b € {2,3,4}, which appears in [4, Theorem 4.2]. They accomplish 
this via the use of a different region 7, than we chose but which achieves the same 
result. In this sense, Theorem 1.4 can be viewed as an extension of [4, Theorem 4.2] 


to b > 20. 


Proof. In addition to the inequalities in (3.9), we note that 


T T 


Pel) = | rotan(/3/(2b + i cinty ey) 


= D(b). 


Unlike (3.9), the above is not a strict inequality. When equality occurs is discussed 
in more detail after this proof. 


Recall the proof of Lemma 3.3. We set 
(€2, €3, €4, €6, d) = (0, i 1, 1, 1.) 


AT 


We define F;(z) as in (3.1), P,(z, y) as in (3.2), and Ry, as in (3.4). In addition to 
D4 = D,4(b), Ds = D3(b), De = De(b), and D = D(b), we set for b € [5,26] the 
values 0 = V(b) = 7/D3(b) and m = m(b) = 1732/(1000(2b — 1)). For b > 26, we set 
0 = V(b) = 7/D(b) and m = m(b) = 1732/(1000(2b + 1)). We note that for each b, 
the number m is a rational number. 

We consider the line y = tan(v)a or equivalently the points x + itan(v)ax in the 


complex plane. Using our choices for ? and m, we obtain for b > 26 that 


tan(¥) =t is a ( ; ( 1732 )) 
ee an | arctan | —————- ] ] = m 
| camo | 1000(26 + 1) 


A calculation shows that tan(v) > m also holds for b € [5,26]. So the line y = ma 


lies strictly below the line y = tan(V)x for x > 0. Applying Lemma 3.3, we know 
that p,(b — ao) = 0 and that p,(x) is continuous. We claim that y = mz lies entirely 
above the region R, for b > 26. We use a Sturm sequence to verify that when b = 27, 
the polynomial P27(7,m(27)x) has no real roots. Using Lemma 3.3 (ii) and (iv), we 
deduce that the line y = m(27)x does not intersect the region Ro7. Using a Sturm 
sequence with P27(x, 1732/2000), we deduce similarly that the line y = 1732/2000 
does not intersect Ro7. Since Ry, is simply a horizontal translation of R27, we see 
that the line y = 1732/2000 lies strictly above the region 7, for all b > 27. The 
lines y = 1732/2000 and y = mz intersect at the point (b+ 1/2, 1732/2000). The 
smaller angle between these two lines is arctan(m/(b)), which is < arctan(m(27)) for 
b > 27. It follows that, since the line y = m(27)z lies strictly above the region Roz, 
we have also that for b > 27, the line y = m/(b)z lies strictly above the region Rp. 
For b € [5,26], we verify similarly with a Sturm sequence that y = m(b)z lies strictly 
above the region Rp. 

We recall the set-up from Section 3.2. We suppose f(x) is reducible and write 
f(x) = g(x)h(x), where both g(x) and A(x) are in Z[x], g(x) F £1, h(x) # +1, and 
both g(x) and h(x) have positive leading coefficients. Furthermore, without loss of 


generality, we suppose that g(b) = +1. In Section 3.2, we showed that either g(x) has 
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a root in common with at least one of 3(a—b), ®4(a—b), and ®g(x—), or g(x) has a 


root 6 € Ry. Since f(x) has non-negative coefficients and the real numbers in Fy are 


positive, we know that 6 ¢ R. So take 6 = 0 + it such that 0 <t < mo < tan(W)o. 
Note the latter implies that we can alternatively write 6 = re”, with r > 0 and 
One ey: 

By Theorem 3.5, we have that f(x) is irreducible if deg f(x) < D4(b). We now turn 
to establishing the statements concerning D3 and Dg. In order to apply Lemma 3.2, 
we observe that 


0<0 <= 5 for b € [5, 26] 
3 


and 


0<0'<d= F< for b> 27. 
3 


For b > 5, we have from (3.9) that 


arg(b+ Cs) < arg(b+ 3) < -_ 
3 


If f(x) is reducible and deg f(x) < D3, we know that f(x) cannot have a root in 
Ry, since otherwise f(x) has a root 8 = re’ € Ry, with r > 0 and 0 < W < a/Ds, 
contradicting Lemma 3.2. Similarly, since arg(b + Gs) < arg(b + ¢3) < m/D3, the 
polynomials ®3(z — b) and ®¢(a — b) cannot be factors of f(x), so f(x) is divisible 
by ®4(a — ). 

For b > 26, we have from (3.9) that arg(b + ¢s) < a/Dg. Thus, we know that 
if f(x) is reducible and the degree of f(x) is < Dg, then from Lemma 3.2 and 
0< WwW <a/D < m/D¢ we obtain that f(z) does not have a root in R, and f(x) is 
not divisible by ®g(x — b), so that f(a) must be divisible by ®3(a — b) or ®4(a — 5). 
Last, for b > 26, we have from Lemma 3.2 and 0 < J’ < a/D that if f(x) is reducible 


and the degree of f(x) is < D, then f(x) cannot have a root in Ry, so that f(x) is 


divisible by @3(a — b), ®4(x — b), or Pg (ax — B). 
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The values D and Dg coincide for finitely many values of b. They first differ when 
b = 64 and they last agree when b = 9355. These facts can be shown using the 


notation of Lemma 3.7, where 


s(b) = 5 


6 


7 500/3(2b4+1) 1 
and thus D(b) = ss ( 1739 — 5) 


To show these values differ when 6 is large, consider the derivative dilg/db. Using 


(3.11), we have for 6 > 2 that 


Ag 


nm Im (b+ C6) 


db |b + Cel? arctan? (Im (b + G6) / Re (b + G)) 


m Im (b+ G6) - 1+2y/1 + (Im (6+ G) / Re (b+ G))” 
Jb + Cel” Im (b+ Cg) / Re (6 + G6) 
~ Tn bt e) (Reb ta) rbr aly’ 
eS bSeal Im (b + C6) 
an 2|b+ 6 ) = Bo a LAY 
Jb + G6? mIim(b+)/ wIm(b+¢.)  2V3 


We see then that Ug(b) is increasing for b > 2 and, by the Mean Value Theorem, 
Ug(b +1) > Ug(b) + 1 for b > 2. Thus, to obtain D(b) > D6(b), it suffices to show 
that 


(OnE 55 
Solving, we see this inequality holds for b > 34088. One checks directly that for all 
integers b € (9355, 34088], we have D(b) > Dg(b). 

Via a generalization of the argument used by Filaseta in [7], we can show that the 


bounds D3, D4, and Dg appearing in Theorem 1.4 are sharp, i.e., they are the largest 


values having the stated properties. 


Theorem 3.8. Letb>5€ Z. For each n € {3,4,6}, let mn, € Z*, Mn > Dp(b) +1. 
Then there are an infinite number of polynomials f,(x) € Zlx| of degree m,, with 


non-negative coefficients for which f,(b) is prime and ®,(ax — b) | fr(x). 
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We note that this result has already been shown for D, in Lemma 3, Lemma 4 
and Theorem 8 of |[7|. The proof below reproduces those arguments with enough 


generalization to show the same result for D3 and Dg. 
Proof. Fix 
ps =(P —b4+1)"?, p= (P41), pe =(P +041)? 


and note that for each n € {3, 4,6} we have 


®,,(x _ b) = (x _ kK) ees = ) where K= pne™, —— pne en. 


Set s, = D,(b) — 1. Since b > 5, we have that s, is non-negative. For j € 
{—1,0,1,...,5n,5n +1}, define 
Kon—j+l _ \on—Jt1 


C= 
: KA 


Observe that the numbers defined by d; = d;(n) = Cs, 41-j,n Satisfy do = Cs,,41n = 9, 


d¢—cy..4° = Land 


diy = Cp pe = (BPA Ca BAC aie 


= (kK =F A)dj41 — KAd; for 0 < j < Sy. 


The above is a recursion relation for d; with characteristic polynomial ®,,(a — b). We 
deduce that c;, € Z for each j € {—1,0,1,...,5n,5n +1} and each n € {3, 4,6}. 


Furthermore, setting 
Sn Sn Snt1 
_ i oe erage 
Un(@) =D) Cjnt? = Yes, jnB 7 = DP dae, 
j=0 j=0 j=l 
we see that the above recursion implies 
ont _ \Snt2 Kent _ Monti 


= — 7Snt2 _ \ 
®, (x2 — b)un(x) = x rs eae =a KA. 


Then since cjn = ps"! sin((Sn — 7 + 1)0n)/sin(6,), we obtain 


Sn+2 gph] sin((Sn + 2)6n) L 8n+2 sin((S, + 1)6n) 
— Pn F vT Py : 
sin(,,) sin(@,) 
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Recalling (3.5), we have 
0 < (Sn +1)On = DnOn < (7/On) + On = 7 < (Sn + 2)On = Dn On + On < 17+ On < 2r. 


We obtain that sin((s, + 2)0,) is negative and sin((s, + 1)6,) is positive, so that the 
coefficient of x and the constant term in ®,,(a — b)u,(x) are positive. 
Taking m, as in the statement of the theorem, observe that we have my, > S, +2. 


To obtain a polynomial of degree m,, we define 
hin (x) = Un (x) (a 897? 4 genom} tt 1) 


so that for an n € {3,4,6}, ®,(x” — b)h, (x) is a polynomial of degree m,, having all 


non-negative coefficients. Then set 
Wy(x) = (20+ 1l)hy(xz) +1. 


Since the coefficients of ®,(x” — b)h,(a) are non-negative and the coefficient for x is 


non-zero, the coefficient of x in ®,(x — b)h,(x) must be at least 1. This implies that 
®,,(a — b)wp(x) = (2b + 1)®,,(x — b)h, (x) + ©, (a — 0) 


has non-negative coefficients, as the least negative coefficient appearing in ®, (ax — b) 


for each n € {3,4,6} is —1 — 2b when n = 6. Furthermore, the values 
hn(b) = ®,(b — b)h,(b) and w,(b) = (2b+ 1)h,(b) +1 


are both positive and relatively prime. Thus every sufficiently large integer is of the 
form Aphy(b) + YnWn(b) for some a@,,Y, € Zt. In particular, for each n € {3,4,6} 
there are infinitely many primes p, such that for some positive integers Q,,Yn, we 


have Qnhn(b) + YnWn(b) = Dn. Set 
fr(@) = On(t — b)(Anbin(@) + YnWn(@)). 


Then f,(a) has degree m,, is reducible, has all non-negative coefficients, and is such 


that f,(b) = pn, a prime. Since there are infinitely many choices for the prime pp, 
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we get that there are infinitely many choices for f,,(x), completing the proof of the 


theorem. 


The bounds in Theorem 3.5 and Theorem 1.4 can be generalized to an arbitrary 
degree if one allows for the possibility of f(a) being divisible by some polynomial 
taken from a finite set of polynomials. To show this result, we first prove the following 


lemma. 


Lemma 3.9. Fixe € (0,1) and G € Zt. Let W = W(e,G) denote the set of 
g(x) € Zia] for which |g(0)| = 1, deg g(x) € [1,G], and |a| > € for each root a € C 
of g(x). Then W is a finite set. 


Proof. Let w(x) € W, and set n = degg. The definition of W implies n < G. Let 
Q1,-..,;Q@ denote the roots of w(x), appearing to their multiplicity. Let a denote 
the leading coefficient of w(x). Since w(x) € W, we obtain |w(0)| = 1 and, for each 


7 € {1,2,...,n}, we have |a;| > ¢. Hence, 


w(0 1 
Ja| |a| é 
Furthermore, for each j € {1,2,...,n}, we deduce 
= }ar| Jaa]: ++ lon 1/a 1 
}a;| = = en-1 = Gi? 
Jar] [aa] ++ faj—1] loj4i]-++len-i| lan] ~ € Jale 


As each coefficient of w(x) after the leading coefficient is +a times a sum of a product 
of distinct a;, we see that the absolute value of each coefficient of w(a) is bounded 
above by 
jal2" max{1/(ale°~*), 1}" < |al2° max{1/(|ale"), 1}° 
Ja|2° G 
aq + lal2 


~ (lale?~*) 


ie Oe oe 
= -(G-NG a eG" 
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The lemma follows as we now see that the degree of w(x) is bounded and the absolute 


value of each coefficient of w(x) is bounded by a function depending only on the fixed 


numbers ¢ and G. 


Theorem 3.10. Let b and G be integers > 2. Then there is a finite set S = S(b,G) 
of polynomials in Z\a| such that, for all f(x) € Zax] with f(x) having non-negative 
integer coefficients, f(b) prime and deg f < G, either f(x) is irreducible or f(x) is 


divisible by some polynomial g(x) € S. 


Proof. It suffices to establish the theorem for b > 2 fixed and G sufficiently large, so 
we take G so that 


bsin(7/G) € (0,1). 


Suppose f(z) € Zl] with f(x) having non-negative integer coefficients, f(b) prime, 
deg f < G, and f(x) reducible. Since f(b) is prime, we have that f(x) = g(x)h(z) 
with g(x) and A(x) in Zz], with deg g(x) > 1, and with g(b) = £1. Since the roots 
of g(x) are roots of f(z) and deg f(x) < G, Lemma 3.2 implies that g(x) cannot 
have any root a = re” with r > 0 and 0 < @ < ma/G. One checks that the line 
passing through the origin with slope tan(7/G) is tangent to the circle centered at 
b of radius bsin(7/G). Taking « = bsin(z/G) € (0,1) and w(x) = g(a +b) in 
Lemma 3.9, we deduce that there are a finite number of possibilities for the factor 


g(x) of f(a). Taking S to be the set of these finitely many possibilities for g(x), the 


theorem follows. 
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APPENDIX A 


COVERINGS FOR BASES OTHER THAN 10 


For bases other than 10, we wish to establish the statement analogous to Theo- 
rem 1.3. To do so, we exhibit coverings of the integers corresponding to the following 


generalization of Lemma 2.2. 


Lemma A.1. Let A,b € Z*, and for a prime p relatively prime to b let cy(p) be 
the multiplicative order of p mod b. For a fixed 6 € {0,...,b—1}, suppose we have 


distinct primes py,..., P14, each relatively prime to b, satisfying the following. 
(i) There exists a covering of the integers 


k=b; (mod @(p;)), 1<i<t. 


(ii) The number A satisfies each of the congruences 


A=-6-b% (mod p,), 1<i<t. 


Then, for each k € Z* U {0}, the number 
A+6-b* 
is divisible by at least one of the primes p; where 1 <i <t. 


The proof for Lemma A.1 is identical to that of Lemma 2.2. In the tables that 
follow, we have at times truncated the presentation of the congruences chosen so as to 
meet formatting guidelines. For example, in Table A.12 we have removed the columns 


with heading ‘row’ to allow for more horizontal space and have also shortened the 
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presentation of each residue class chosen. Whenever possible, each covering progresses 
roughly for some m from k = 0 (mod m) through k = m-—1 (mod m). However, we 
have at times rearranged the rows in the interest of placing any lengthy primes in a 
single column. While this presentation could make the original logic of the covering 


difficult to follow, it serves to keep white space from over-abounding in the tables. 


A.1 COVERING SYSTEMS FOR b = 2 


For b = 2, our system of congruences results in a modulus M on the order of 10!%, or 
2°" In Table A.1 we exhibit a covering corresponding to Lemma A.1 (i) which allows 
a leading 0 in binary to be changed to a 1. 


Table A.1 Covering for A+1- 2* 


row congruence prime p; 
1 k =0 (mod 2) 3 
2 k =1 (mod 4) 5 
3 k =3 (mod 8) Alb 
4 | k=7 (mod 16) | 257 
5 | k=15 (mod 32) | 65537 
6 | k=81 (mod 64) 641 
7 | k=63 (mod 64) | 6700417 


A.2 COVERING SYSTEMS FOR b=3 


For b = 3, our system of congruences results in a modulus M on the order of 10%, 


or 371". The smallest widely digitally delicate probable prime in base 3 that we have 


found using such a system is 
5120877327448 1310348328072804273775352870316957766600352463048 136275185— 


6722041264394581318449651152166653. 


We take A = 1 (mod 2) so that A+ 1- 3* is divisible by 2. Table A.2 exhibits a 


covering which allows for a leading 0 to be changed to a 2 in base 3. 
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Table A.2 Covering for A+ 2-3* 


congruence prime p; congruence prime p; 

k =0 (mod 3) 13 k =17 (mod 72) 282429005041 

k =1 (mod 6) 7 k = 53 (mod 144) 1418632417 

k =2 (mod 9) 757 k = 125 (mod 144) 56227703611393 

k =8 (mod 18) 19 k = 23 (mod 216) 2161 
k=14 (mod 18) | _37 k = 95 (mod 216) 15121 

k =0 (mod 4) 5 k = 167 (mod 216) 10512289 
k = 10 (mod 12) 73 k = 71 (mod 108) | 150094634909578633 
k =5 (mod 36) | 530713 k = 143 (mod 216) 16569793 

k =3 (mod 8) Al k = 215 (mod 216) 3958044610033 


A.3 COVERING SYSTEMS FOR }6 = 4 


For b = 4, our system of congruences results in a modulus M on the order of 10°°, or 
4°°. We take A = 2 (mod 3) so that A+1- 4* is divisible by 3. Tables A.3 and A.4 
exhibit covers which allow for a leading 0 to be changed to a 2 or a 3 in base 4. 


Table A.3 Covering for A+ 2- 4* 


row congruence prime p; 
1 k =0 (mod 2) 5 
2 k =1 (mod 4) 17 
3 k = 3 (mod 8) 251 
4 k =7 (mod 16) 65537 
5 | k=15 (mod 32) 641 
6 | k=83l (mod 32) | 6700417 


Table A.4 Covering for A+ 3- 4* 


row congruence prime p; 

1 k =0 (mod 3) 7 

2 k =1 (mod 6) 13 

3 | k=4 (mod 12) 241 
4 | k=10 (mod 24) 97 
5 | k = 22 (mod 24) 673 
6 k =2 (mod 9) 19 
i k =5 (mod 9) 73 
8 | k=8 (mod 18) at 
9 | k=17 (mod 18) 109 
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A.4 COVERING SYSTEMS FOR b=5 


For b = 5, our system of congruences results in a modulus M on the order of 10!%, 


or 52. We take A = 1 (mod 2) so that A+1-5* and A+3-5* are divisible by 


2. We exhibit coverings in the tables below corresponding to the other possible digit 


increases. 


Table A.5 Covering for A+ 4-5* 


congruence prime p; congruence prime p; 
k = 0 (mod 5) 11 k = 23 (mod 150) 118801 
k =1 (mod 5) 71 k = 53 (mod 150) 20775901 
k = 3 (mod 6) 7 k = 83 (mod 150) 24665701 
k =2 (mod 15) 181 k = 113 (mod 150) 149439601 
k=7 (mod 15) | 1741 k = 143 (mod 300) 14401 
k=4 (mod 10) | 521 k = 293 (mod 300) | 299541552154912341601 
k =8 (mod 20) Al k = 29 (mod 90) 60081451169922001 
k=18 (mod 20) | 9161 k =59 (mod 180) 20478961 
k=12 (mod 30)| 61 k = 149 (mod 180) 6794091374761 
k =13 (mod 30) | 7621 k = 89 (mod 180) 25535754811081 
k = 19 (mod 60) 2281 k = 179 (mod 360) 8641 
k =49 (mod 60) | 69560521 | [| & =359 (mod 360) 140641 
Table A.6 Covering for A +2 -5* 
congruence prime p; congruence prime p; 
k=0(mod3) | 31 k = 11 (mod 24) 390001 
k =0 (mod 2) 3 k = 23 (mod 48) 152587500001 
k =1 (mod 6) 7 k = 47 (mod 96) 97 
k =5 (mod 12) 601 k = 95 (mod 96) | 240031591394168814433 


A.5 COVERING SYSTEMS FOR b=6 


For b = 6, our system of congruences results in a modulus M on the order of 1 


483 
aaa 


T 


6°70. We take A = 2 (mod 5) so that A+3-6* is divisible by 5. We exhibit coverings 


in the tables below corresponding to increasing a digit in base 6 by the other possible 


increases of 1, 2, 4, or 5. 
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Table A.7 Covering for A+ 1- 6* 


row congruence prime p; 
rice aoe) a row congruence prime p; 
= 6 | k=3 (mod 18) | 46441 
2 | k=0 (mod 9) 19 — 
- 7 | k =12 (mod 36) 73 
3 | k =6 (mod 9) 2467 = 
— 8 | k =30 (mod 36) 541 
4 | k=1 (mod 12) 13 9 =O (aniod'S) B 
5 |k=7 (mod 12)|__97 ao 
Table A.8 Covering for A+ 2- 6* 
row | congruence | prime p; 
Se = row congruence prime p; 
— 4 | k=7 (mod 16) 17 
see) ee 5 | k=15 (mod 16) | 98801 
3 | k=3 (mod 8) 1297 — 
Table A.9 Covering for A+ 4- 6* 
congruence prime p; ieee ese ae 
(ied 5) 311 58 (mod 60) 181 
Tcaed10) = 4 (mod 25) 18198701 
a ae) 7 9 (mod 25) 40185601 
14 (mod 50) 3655688315536801 
2 (mod 15) 1171 
Ted = 39 (mod 100) 343801 
89 (mod 100) 22243201 
12 (mod 45) 2161 
19 (mod 75) 601 
27 (mod 45) | 112771 
42 (mod 45) | 19353635731 znd) eal 
3 (mod 20) Al 69 (mod 75) 271041511600591342728451 
ae Aan 24 (mod 125) 9536585501 
13 (mod 40) 7 49 (mod 125) 117811792772681609501 
33 (mod 40) | 68754507401 74 (mod 125) | 105875321588567599765751 
99 (mod 125) | 1098445767808750903973251 
18 (mod 60) 61 
eV aaloDen Sel 124 (mod 250) 251 
249 (mod 250) 751 


We use each of the primes 7, 13, and 31 twice over Tables A.7, A.8 and A.10. 


One checks that our choices are equivalent to A= 5 (mod 7), A=7 (mod 13), and 


A=6 (mod 31). For Table A.10 we set 
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Pir = 2245225770735455723534882978547 1057921. 


Table A.10 Covering for A+ 5 - 6” 


congruence prime p; congruence prime p; 
1 (mod 2) 7 54 (mod 108) 591841 
2 (mod 6) 31 90 (mod 108) 171467713 
4 (mod 12) 13 24 (mod 108) 932461936453 
10 (mod 24) 1678321 60 (mod 216) 133 
22 (mod 48) 5953 168 (mod 216) 115963921 
46 (mod 48) 473896897 96 (mod 216) 8781208996949976153601 
0 (mod 36) 55117 204 (mod 216) | 241282001155985351966017 
6 (mod 72) 577 30 (mod 180) 9001 
42 (mod 72) 3313 66 (mod 180) 211501 
12 (mod 72) | 2478750186961 102 (mod 180) 2106930961 
48 (mod 144) Pi 138 (mod 180) 5997780112726834061 
120 (mod 288) 115777 174 (mod 540) 4861 
264 (mod 288) 31057921 354 (mod 540) 39326041 
18 (mod 108) 109 534 (mod 540) 51353541541 


A.6 COVERING SYSTEMS FOR b = 7 


For b = 7, our system of congruences results in a modulus M on the order of 10!, 
or 7841. We take A = 1 (mod 2) so that A+ 1-7*, A+3-7*, and A+5-7* are 
divisible by 2. We also take A = 1 (mod 3) so that A +2-7* is divisible by 3. We 
exhibit coverings below corresponding to the other possible digit increases in base 7. 


For Table A.12, we set 


Piz = 217648180992721729506406538251, pis = 2182816753758823696751, 


Pig = 9903546678356844440204179342119473011887617239565323800435501, 
pag = 271796439196451766191391111220539009912588742400633481231401, 


pag = 353207957997402826578121, — pso = 11193638545408118035815864041. 


We obtain the factorization 
®499(7) = 491 - 12251 - 20632921 - 648128907712931 - pso - Caso 
and again use the notation pn1, Pn to refer to two distinct primes dividing C’,. 
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Table A.11 Covering for A+ 4-7* 


TOW congruence prime Di Ou: Ouse ne Fs 
; : - : a : = 7 | k=17 (mod 24) 193 
3 k = 1 Gc 6) 43 8 k=11 (mod 24) 409 
=o k = 23 (mod 48) | 33232924804801 
4 k =2 (mod 12) 13 = 
5 |k=10 (mod 12)|_ 181 10_| k = 47 (mod 96) 97 
= = 104 
6 | k=5 (mod 24) 73 11 | k =95 (mod 96) 04837857 
Table A.12 Covering for A +6-7* 
: congruence prime p; 
congruence prime p; 
0 (mod 5) 801 57 (mod 60) 555915824341 
I (mod 10) Tl 8 (mod 25) 31280679788951 
6 aed 10) 19 13 (mod 50) 79787519018560501 
5 ca 15) 3] 143 (mod 150) 6005492312551 
7 ae 15) 159871 223 (mod 250) 1173001 
1D a ~ 30) | 6568801 123 (mod 250) 4833574921448501 
— 248 (mod 250) 31818863467130251 
at mot = 4 (mod 35) 2127431041 
3 (mod 25) 2551 
9 (mod 35) 77192844961 
38 (mod 100) 101 
19 (mod 70) 12128131 
88 (mod 100) | 13001 
18 (mod 75) | 29251 54 (mod 70) 603926681 
43 (mod 75) 24 (mod 140) 102314938321 
68 ned 150) Hi 94 (mod 140) 18690488255321 
33 aed 125) al 129 (mod 210) | 338640865331157691 
7B fied 125) 199 (mod 210) | 450798894542150330401 
~“ oe 64 (mod 175) | 1237578612719152201 
73 (mod 125) 016 
98 (mod 250) | 751 99 (mod 175) Dag 
14 (mod 70) A 134 (mod 175) | 1094471317606762277701 
19 aed 70) i 344 (mod 350) 8230203760252601 
59 nied 210) a1 34 (mod 245) 4896864001 
99 aed 175) 701 69 (mod 245) 9152363081 
169 ee ; d 350) | 3851 104 (mod 245) 11492120512321 
139 eed 245) 209 (mod 490) 20632921 
“ es 454 (mod 490) 648128907712931 
174 (mod 490) | 491 
419 (mod 490) | 12251 244 (mod 490) Ds0 
“ 489 (mod 490) P4901 
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A.7 COVERING SYSTEMS FOR b=8 


For b = 8, our system of congruences results in a modulus M on the order of 101778, 
or 8'%!4. We take A = 4 (mod 7) so that A+ 3- 8* is divisible by 7. We exhibit 
coverings in the tables below corresponding to the other possible digit increases in 
base 8. 

We make repeated use of the primes 3, 5, and 17 in Tables A.13, A.14, A.15 and 


A.16. One checks that our choices are equivalent to A= 1 (mod 3), A=4 (mod 5), 


and A= 4 (mod 17). 


Table A.13 Covering for A+1-8* 


row congruence prime p; row | congruence | prime p; 
1 k =1 (mod 2) 3 3 | k=2 (mod 8) 17 
k =0 (mod 4) 5 4 | k =6 (mod 8) 241 
Table A.14 Covering for A +2-8* 
row congruence prime p; 
1 k =0 (mod 2) 3 
2 k =1 (mod 4) 5 
3 k =3 (mod 4) 13 
Table A.15 Covering for A+ 4-8* 
row | congruence | prime p; 
ieee) 7 row congruence prime p; 
= 4 | k=0 (mod 16) 97 
2 WSemod a 5 |k=8 (mod 16) | 257 
3 |k=4(mod8)| 17 aa 
Table A.16 Covering for A+5-8* 
row congruence prime p; 
i FE ood) 7 row congruence prime p; 
= 5 | k=7 (mod 10) | 331 
2 | k=0 (mod 5) 31 — 
— 6 | k=9 (mod 20) A] 
3 | k=1 (mod 5) 151 7 |k=19 (mod 20) 61 
4 |k=3(mod10)| 11 = 


For Table A.20, we list partial factorizations of ®,(8) in Table A.18. 
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Table A.17 Covering for A+ 6 - 8* 


row congruence prime p; ; 
= A 3) 73 TOW congruence prime p; 
: : = : a“ 16) 7 7 |k=14 (mod 36) | 246241 
= 8 | k =32 (mod 36) | 279073 
3 | k =4 (mod 12) 37 — 
= 9 | k=8 (mod 27) | 2593 
4 | k =10 (mod 12) 109 = 
10 |k=17 (mod 27) | 71119 
I cto) 20202 1 | k =26 (mod 27) | 97685839 
6 | k=5 (mod 18) | 87211 = 


Table A.18 Partial factorizations of 


®,,(8) for large n 
n | Partial factorization of ®,,(8) 
245 | 1471-41161 - 4163041 - Co45 
343 | 6073159 - 1428389887 - C343 
686 8233 - 2513690593 - Cege 
1029 271657 - Co29 


Table A.19 First part of covering for 


A+7-8* 
row congruence prime p; 
1 k =0 (mod 7) 127 
2 k =1 (mod 7) 337 
3 | k =2 (mod 14) 43 
4 | k=9 (mod 14) 5419 
5 | k=3 (mod 21) 92737 
6 | k=10 (mod 21) 649657 
7 | k=17 (mod 42) | 77158673929 
8 | k =38 (mod 84) | 40388473189 
9 | k =80 (mod 84) | 118750098349 
10 | k =4 (mod 35) 71 
it |k=11 (mod 35) 29191 
12 | k=18 (mod 35) 106681 
13 | k = 25 (mod 35) 122921 
14 | k = 32 (mod 35) 152041 
15 | k=5 (mod 70) 214 
16 | k=12 (mod 70) 281 
17 | k =19 (mod 70) 86171 
18 | k = 26 (mod 70) 664441 
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Table A.20 Second part of covering for A+ 7 - 8* 


row congruence prime p; row congruence prime p; 

19 | k= 33 (mod 70) 1564921 42 | k =69 (mod 98) 4363953127297 

20 | k =40 (mod 140) 421 43 | k =27 (mod 98) 2603288584539209385 1 

21 | k=110 (mod 140) 7416361 44 k = 76 (mod 196) 197 

22 | k =47 (mod 140) 146919792181 A5 k = 174 (mod 196) 540961 

23 | k = 117 (mod 140) 1041815865690181 46 k = 34 (mod 196) 19707683773 

24 | k=54 (mod 210) 1765891 AT k = 83 (mod 196) 4981857697937 

25 | k = 124 (mod 210) 1124770259967650548- 48 k = 132 (mod 196) 40544859693521 152369 
144713799166434869 1 AQ k = 181 (mod 196) | 17059410504738323992180849 

26 | k = 194 (mod 420) 2521 50 | k&=41 (mod 245) 1471 

27 | k = 404 (mod 420) 1711081 51 | &=90 (mod 245) 41161 

28 | k= 61 (mod 350) 1051 52 | k= 139 (mod 245) 4163041 

29 | k = 131 (mod 350) 110251 53 k = 188 (mod 245) P245,1 

30 | k = 201 (mod 350) 3205651 54 | k = 237 (mod 245) Pais 

31 | k = 271 (mod 350) 247772800801 50 k = 48 (mod 343) 6073159 

32 | k = 341 (mod 350) 347833278451 56 k = 97 (mod 343) 1428389887 

33 | k =68 (mod 350) 7223591273619001 oT k = 146 (mod 343) P343,1 

34 | k = 138 (mod 350) 34010032331525251 58 k = 195 (mod 343) 2343,2 

35 | k = 208 (mod 350) 129266711542799251 59 k = 244 (mod 686) 8233 

36 | k = 278 (mod 350) 2310141222312973778401 60 k = 587 (mod 686) 2513690593 

37 | k = 348 (mod 700) 701 61 | k= 293 (mod 686) P686,1 

38 | k = 698 (mod 700) 6301 62 | k = 636 (mod 686) P686,2 

39 k =6 (mod 49) 4432676798593 63 | k = 342 (mod 1029) 271657 

AO k = 13 (mod 49) | 2741672362528725535068727 64 | k = 685 (mod 1029) P1029,1 

41 | k =20 (mod 98) 748819 65 | k = 1028 (mod 1029) P1029,2 


A.8 COVERING SYSTEMS FOR b = 9 


For b = 9, our system of congruences results in a modulus M on the order of 101”, 
or 9188. We take A = 1 (mod 2) so that A+ (27 +1)-9* is divisible by 2. We exhibit 


coverings below corresponding to increasing a digit in base 9 by an even amount. 


We use each of the primes 5 and 17 twice over Tables A.22, A.23, and A.24. One 


Table A.21 Covering for A+ 6-9* 


row congruence prime p; 
1 k =0 (mod 3) ec 
2 k =1 (mod 3) 13 
3 k =2 (mod 6) 73 
4 | k=5 (mod 12) 6481 
5 | k=11 (mod 24) 97 
6 | k= 23 (mod 24) 577 


checks that our choices are equivalent to A= 3 (mod 5) and A= 4 (mod 17). 


Table A.22 Covering for A +2-9* 


row congruence prime p; 
i k =0 (mod 2) 5 
2 k =1 (mod 4) Al 
3 | k=3 (mod 8) 17 
4 | k=7(mod8) | 193 
Table A.23 Covering for A+ 4-9* 
anid eee ere puma? row congruence prime p; 
pil 2 See ©) a 9 |k=12 (mod 15)| 4561 
2 k = 0 (mod 5) 11 = 
3 eA nod 5) i 10 | k =4 (mod 25) 151 
7 = (mod 10) Tisi 11 | k=9 (mod 25) 8951 
= 12 |k=14 (mod 25) | 391151 
5 | k=18 (mod 20) | 42521761 = 
= 13 | k=19 (mod 25) | 22996651 
6 | k=8 (mod 40) | 14401 = 
7 = 2 mod 1s) 3] 14 | k = 24 (mod 50) 101 
mere 15 | k =49 (mod 50) | 394201 
8 | k=7 (mod 15) | 271 
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Table A.24 Covering for A+8-9* 


row congruence prime p; 
1 k =1 (mod 2) 5 
2 k = 0 (mod 7) 547 
3 | k=1 (mod 7) 1093 
4 k =2 (mod 14) 29 
5 | k=4 (mod 14) 16493 
6 k =6 (mod 28) 430697 
7 k = 20 (mod 28) 647753 
8 k = 10 (mod 42) 2857 
9 | k=24 (mod 42) 109688713 
10 | k =88 (mod 84) 337 
11 | k =80 (mod 84) 673 
12 k = 12 (mod 70) 28596961 
13 k = 26 (mod 70) 32839661 
14 k = 40 (mod 70) 94373861 
15 | k= 54 (mod 140) 281 
16 | k = 124 (mod 140) 18481 
17 | k =68 (mod 140) 369879560116990841 
18 | k = 138 (mod 140) | 3353336738929580410561 


A.9 COVERING SYSTEMS FoR b= 11 


For 6 = 11, our system of congruences results in a modulus M on the order of 1 
or 11°!. We take A = 1 (mod 2) so that A+ (2j +1)-11* is divisible by 2. We also 
take A = 3 (mod 5) so that A +2-11* is divisible by 5. We exhibit coverings below 
corresponding to increasing a digit in base 11 by other possible amounts. The prime 


3 is used multiple times and in each case corresponds to A = 2 (mod 3). 


Table A.25 Covering for A+4-11* 


row congruence prime p; 
1 k =0 (mod 2) 3 
2 k =1 (mod 4) 61 
3 k =3 (mod 8) 7321 
4 | k=7 (mod 16) 17 
5 | k=15 (mod 16) | 6304673 


69 


Table A.26 Covering for A+ 10-11* 


: congruence prime p; 
congruence | prime p; 
25 (mod 56) 77001139434480073 
0 (mod 2) 3 
0 (mod 7) B 39 (mod 112) aot 
Pisa) ss 95 (mod 112) 394129 
3 (mod 14) 1623031 53 (mod 112) 236352238647181441 
5 (mod 28) 29 109 (mod 112) | 3090443962383595123379137 
13 (mod 42) 3421169496361 
19 (mod 28) 1933 
Gi(nad 28) | Bopa7aTg || | med Be) ae 
69 (mod 84) 520799717831587692709 
23 (mod 56) 113 
oe ee 1 (mod 126) 3304981 
iad $6) O52 83 (mod 126) 468843103 
125 (mod 126) 71596275661 
Table A.27 Covering for A+ 8- 11* 
congruence prime p; a prime p; 
I (mod 2) 3 4 (mod 40) 1120648576818041 
0 (mod 5) 3221 6 (mod 90) 86306335830799838011 
2 (mod 10) | 13421 8 (mod 50) 14436295738510501 
4 (mod 20) | 212601841 7s (mod 100) 2248313994601 
14 (mod 40) Al 38 (mod 100) | 1993099906542710819727884501 
6 (mod 30) 31 88 (mod 200) 401 
16 (mod 30) 7537711 188 (mod 200) 2254601 
26 (mod 90) 181 48 (mod 150) 751 
56 (mod 90) 631 98 (mod 150) 299401 
8 (mod 50) | 46601 148 (mod 150) 20128755258227254- 
28 (mod 100) 101 47651340174136551 
Table A.28 Covering for A+6-11* 
row | congruence | prime p; 
i. 1F=0. Giedg) = row congruence prime p; 
— 4 | k=4 (mod 12) 13 
2s |= Ni ods), 5 |k=10 (mod 12)| 1117 
3 |k=2(mod6)| 37 = : 


A.10 COVERING SYSTEMS FOR }) = 31 


For b = 31, we did not determine the size of the modulus MW resulting from our system 


of congruences. We take A = 1 (mod 2) so that A+ (2j+1)-31* is divisible by 2. We 
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also take A = 2 (mod 3) so that A+(3j+1)-31* is divisible by 3, and we furthermore 
set A = 1 (mod 5) so that A+ (5j +4) -31* is divisible by 5. We exhibit coverings 
below corresponding to increasing a digit in base 31 by other possible amounts. The 
primes 7, 11, 13, 17, 19, 37 are used multiple times below and each use corresponds 


to one of the following congruences. 


A=5 (mod7), A=5 (mod1l), A=5 (mod 13), 
A=12 (mod17), A=11 (mod19), A=11 (mod 37). 
Table A.29 Covering for A+ 26 -31* 

row congruence prime p; 

1 k = 3 (mod 6) 7 

2 k =1 (mod 6) 19 

3 k =0 (mod 4) 37 

4 k =11 (mod 16) 17 

5 k =5 (mod 18) 577 

6 | k=11 (mod 18) 1538083 

ic k =17 (mod 36) 1536553 

8 | k=35 (mod 36) 512616735577 

9 | k=2 (mod 12) 922561 

10 k =6 (mod 24) 852890113921 

11 | k = 22 (mod 72) 73 

12 | k =46 (mod 72) 4683817 

13 k = 70 (mod 72) 1814503763676 130449408979921 
14 | k= 18 (mod 120) 35401 

15 | &=42 (mod 120) 1546081 

16 | k= 66 (mod 120) | 9667783133425605119410155998541192601 
17 | k=90 (mod 240) DAI 

18 | & = 210 (mod 240) 32360641 

19 | k=114 (mod 240) 362634922081 

20 | k = 234 (mod 240) 69916284426778163281 

20 | k=10 (mod 48) 97 

20 k = 34 (mod 48) 7499207440683838894753 


For the tables that follow, we list partial factorizations of ®,(31) in Table A.31 


and Table A.32. In addition to the values displayed there, we use the notation 
176(31) = Pize and ®420(31) = Puzo 
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Table A.30 Covering for A+ 30-31* 


row congruence prime p; 
it k =1 (mod 8) 331 
2 k = 0 (mod 6) i 
3 k = 3 (mod 6) 19 
4 k =2 (mod 9) 3637 
5 | k=5 (mod 9) 81343 
6 k = 8 (mod 27) 1836205027201 
7 | k=17 (mod 27) | 126901881805771 
8 | k = 26 (mod 81) 2593 
9 | k=53 (mod 81) 13933 
10 | k = 80 (mod 81) 477739 


for values which were determined to be prime. We also use 
bo9(31) = Coos, P325(31) = C325, Baa5(31) = Cazs, ©520(31) = C20, 


$588(31) = Cg, D1975(31) a Cha75; ®1530(31) = C1530; and © 1540(31) = C1540 


for values which we determined to have at least 2 prime factors but for which we 
could not find any prime divisors. Table A.33 displays complete factorizations with 
notation used for primes where were too lengthy to display. 


For Tables A.35 and A.36 we set 
pq = 27262005427411199665751772615574851310841886637227545549 72281 , 


p31 = 660479435423419861325502790686400761533380000744579973299352935629 — 
42426481, pgg = 224721202412918666334576819250523191369, 
p39 = 348673230951446260771883379412341320615347347601103310801195979086— 
277, ps3 = 1129363636895809892086303692627113871721, 


Pog = 29510535204545262157687088665468 191183896988343413667225871. 
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Table A.31 Partial factorizations of ®,,(31) for large n 


n Partial factorization of ®,,(31) n Partial factorization of ®,,(31) 
91 2549 - 1661479 - Co, 468 937 - 1015561 - Cuge 
117 2534689 - C1417 480 | 20641 - 330991018976905188481 - C49 
125 251 - 129001 - 12181751 - Chas 490 1471 - 32341 - 10996091 - C499 
135 1344742561 - Cy35 095 89278561 - 108837401 - 42498121331 
147 883 + Chaz -3933848803218900604400041 - Cros 
160 380641 - 1176641 - 8084410241 - Che 612 15913 - 34273 - Cero 

195 35956831 - Cig5 650 391301 - 45514951 - 1431069251. 
238 11781795397 - 2744226674742050863 - Co38 680 3576121 - Cégo 
245 491 - 55848731 + Coas 765 1531 - 6121 - 1074061 - 8197741 - C65 
209 934831 - Coss 770 2868251 - Cr79 
264 183516169 - Cogs 780 17329261 - 45796141 - Crg9 
306 | 307 - 78198146102753533 - 19220209997787857 - C'396 850 40801 - Ces 
320 259201 - 59826881 - C329 882 154351 - Cggo 
330 2995081 - C339 910 131041 - Coi0 
340 1361 - 3061 - 30941 - 2964461 - 5153381 - C349 1040 2081 - 2057345681 - Cio49 
301 132679 - 445069 - 175071781 - 1467780211 - C35) 1071 39307843 + C1971 
302 353 + C352 1190 2381 + Ci190 
307 637603 - 2112727 - C357 1224 29841121 - 209530441 - 330188689 
385 2311 - Cg -1133681368703716513 
390 188431581362701 - 9708046814116951 - C99 -61410533730626180017 - Cy204 
392 29401 - 946681 - C92 1248 4993 - 2215201 - Cyo48 
416 1404631073 - C16 1300 4479139601 - C1300 
455 6041004971 - C455 1365 2731 - 2435161 - C365 
459 1131567193 - C59 1470 282506071 - C4709 


Table A.32 More partial factorizations of ®,,(31) for large n 


n Factorization of ®,,(31) 
1560 41461681 - Cys5¢60 
1820 14990810381 - Cig20 
1872 282419281 - Cig72 
2142 | 53551 - 8281014336631 - 17562188355029137 - C142 
2275 50051 - 54601 - Coa75 
2448 637045489 + Co448 
2295 4591 - 9181 - 59671 - C295 
2550 2551 - Co559 
3213 7812268926269851 - C3213 
3744 5814433 - C3744 
4284 5714857 - C284 
4590 2519911 - C4590 


Table A.33 Complete factorizations of ®,(31) for large 


n 
n Factorization of ®,,(31) 
64 A801 - Pea 
85 | 108971 - 391206807721 - 21736504684553261 - Pgs 
119 239 - 2857 - 78541 - Pig 
153 17137 - 42834601810502407 - Pys3 
165 1321 - Pres 
170 732122807601 - Pio 
180 8728381 - 14398921 - Pigo 
196 197 - Pigg 
200 601 - 7137001 - Poo 
260 1301 - 22765081 - Pogo 
312 | 313 - 4292809 - 185272593 - 115436220433 - Ps12 
936 7489 - Po36 


Table A.34 Covering for A+8-31* 


row | congruence | prime p; row congruence prime p; 
1 | k=2 (mod 6) 7 7 | k =83 (mod 16) 17 
2 | k=0 (mod 6) 19 8 k =7 (mod 16) | 25085030513 
3 | k=0 (mod 4) 13 9 | k=11 (mod 32) 1889 
4 |k=1(mod4)| 37 10 | k =27 (mod 32) | 1347329 
5 | k=2(mod8)| 409 11 | &=15 (mod 32) | 6139297 
6 |k=6 (mod 8)| 1129 12 |k=31 (mod 32) | 23277313 
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Table A.35 First part of covering for A +2-31* 


row congruence prime p; row congruence prime p; 

1 k =0 (mod 6) 7 24 k = 33 (mod 91) Poi 

2 k = 15 (mod 16) 17 25 k = 46 (mod 91) Po12 

3 k =3 (mod 5) 11 26 | k=59 (mod 182) 547 

4 | k=0 (mod 13) 42407 27 | k =150 (mod 182) 1093 

5 k =1 (mod 18) 2426789 28 k = 72 (mod 182) 647001 

6 k =2 (mod 18) 7908811 29 | k =163 (mod 182) 14407667 
7 k =3 (mod 26) 17863 30 k = 85 (mod 182) 669665808855863 
8 | k=16 (mod 26) 4276694944587 31 | k=176 (mod 182) Pai 

9 | k=4 (mod 39) 79 32 | k=8 (mod 117) 2534689 
10 | k=17 (mod 39) 13807 33 | k=21 (mod 117) pa 

11 | k=80 (mod 39) 39703 34 | k = 34 (mod 117) Pure 

12 | k=5 (mod 52) 53 35 | k=47 (mod 234) 4447 

13 k = 18 (mod 52) 116337521 36 | k = 164 (mod 234) 7776289 
14 k = 31 (mod 52) 76037563733 37 k = 60 (mod 234) 21931507729 
15 k = 44 (mod 52) 101686136508893 38 | k =177 (mod 234) | 903087677909176579 
16 k =6 (mod 39) 175500339130677572941801 39 k = 73 (mod 234) 139 

17 | &=19 (mod 78) 157 40 | k =190 (mod 468) 937 

18 k = 58 (mod 78) 238213 41 | k = 424 (mod 468) 1015561 
19 k = 32 (mod 78) | 17123370267331917425544180721 42 k = 86 (mod 468) P468,1 

20 | k=71 (mod 156) 141336778441 43 | k = 203 (mod 468) P468,2 

21 | k=149 (mod 156) Pat 44 | k = 320 (mod 936) 7489 

99 k =7 (mod 91) 2549 45 | k = 788 (mod 936) P36 

23 k = 20 (mod 91) 1661479 46 | k = 437 (mod 1872) 282419281 
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Table A.36 Second part of covering for A + 2-31* 


row congruence prime p; row congruence prime p; 

48 | k = 1373 (mod 1872) |__pis7a2 70 | k= 24 (mod 65) 5979236519649901 
49 | k = 1841 (mod 3744) 5814433 71 k = 37 (mod 130) 131 

50. | k = 3713 (mod 3744) P3744, 72 | k=102 (mod 130) 521 

51 k =99 (mod 351) 132679 73 k = 50 (mod 325) 325.1 

52 | k=216 (mod 351) | 445069 74 | k=115 (mod 325) D395,2 

53 k = 333 (mod 351) 175071781 75 k = 180 (mod 650) 391301 

04 k = 112 (mod 351) 1467780211 76 k = 505 (mod 650) 45514951 

55 | k= 229 (mod 351) D351 77 | k= 245 (mod 650) 1431069251 
56 | k = 346 (mod 351) P351,2 78 | k=570 (mod 650) Pe50,1 

Bt k = 10 (mod 65) 911 79 | k=310 (mod 650) P650,2 

08 k = 36 (mod 195) 309096831 80 | k =635 (mod 1300) 4479139601 
59 k= 101 (mod 195) P195,1 81 k = 1285 (mod 1300) P1300,1 

60 | k = 166 (mod 195) P1952 82 | k=12 (mod 65) 31035996941 
61 k = 49 (mod 260) 1301 83 k = 25 (mod 65) ps3 

62 k = 114 (mod 260) 22765081 84 k =9 (mod 104) 305688893141113 
63 k = 179 (mod 260) P69 85 k = 22 (mod 104) 5603212901768856193 
64 | k= 244 (mod 780) | 17329261 86 | k= 35 (mod 104) Dee 

65 | k=504 (mod 780) | 45796141 87 | k =48 (mod 208) pie 

66 k = 764 (mod 780) P780,1 88 k = 152 (mod 208) P208,2 

67 k = 62 (mod 130) 197271101 89 k = 61 (mod 312) 313 

68 | k= 127 (mod 130) Dos 90 | k= 165 (mod 312) 4292809 
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Table A.37 Third part of covering for A+ 2-31* 


row congruence prime p; row congruence prime p; 
91 k = 269 (mod 312) 135272593 113 k = 181 (mod 455) P455,2 

92 k =74 (mod 312) 115436220433 114 k = 246 (mod 910) 131041 
93 k=178 (mod 312) P312 115 k= 701 (mod 910) P910,1 

94 k, = 282 (mod 624) 1249 116 | k =311 (mod 1365) 2731 

95 k = 594 (mod 624) 1873 117 | k = 766 (mod 1365) 2435161 
96 k = 87 (mod 416) 1404631073 118 | k = 1221 (mod 1365) P1365,1 

97 k=191 (mod 416) P416,1 119 k = 376 (mod 910) P910,2 

98 k = 295 (mod 416) P416,2 120 | k = 831 (mod 1820) 14990810381 
99 k = 399 (mod 1248) 4993 121 | k = 1741 (mod 1820) P1820,1 
100 | k = 815 (mod 1248) 2215201 122 | k = 441 (mod 2275) 50051 
101 | k = 1231 (mod 1248) P1248,1 123 | k = 896 (mod 2275) 54601 
102 | k=100 (mod 520) P520,1 124 | k = 1351 (mod 2275) P2275,1 
103 k = 204 (mod 520) P520,2 125 | k = 1806 (mod 2275) P2275,2 
104 | k = 308 (mod 1040) 2081 126 | k = 2261 (mod 4550) 22751 
105 | k = 828 (mod 1040) 2057345681 127 | k = 4536 (mod 4550) 4770770551 
106 | & = 412 (mod 1040) P1040,1 128 | k=64 (mod 390) 188431581362701 
107 | k = 932 (mod 1040) P1040,2 129 k = 129 (mod 390) 9708046814116951 
108 | k = 516 (mod 1560) 41461681 130 k = 194 (mod 390) P390,1 
109 | k = 1036 (mod 1560) P1560,1 131 k = 259 (mod 390) P390,2 
110 | k = 1556 (mod 1560) P1560,2 132 | k = 389 (mod 1170) 2341 
111 k =51 (mod 455) 6041004971 133 | k =779 (mod 1170) 24571 
112 k = 116 (mod 455) P455,1 134 | k = 1169 (mod 1170) 111127771 


For Table A.38 we set 


Pir = 28188789169957630949493 1758688962111082311886906752520001. 


Table A.38 Covering for A+ 6-31* 


congruence prime p; - 
7 (mod 5) A congruence prime p; 
98 (mod 200) Pyo0 
1 (mod 5) 17351 
2 (mod 10) rT 198 (mod 400) 401 
398 (mod 400) 16001 
7 (mod 10) 21821 
4 (mod 25) 101 
3. (ned. 20) a 9 (mod 25) 4951 
13 (mod 20) 4707206941 
5 God 0) = 14 (mod 25) 13277801 
19 (mod 25) | 20235942281002951 
18 (mod 50) 1726301 os = 
28 (mod 50) | 74770514303869505101 
seid 100) — 49 (mod 125) 129001 
= 74 (mod 125) 12181751 
48 (mod 100) Pu 
99 (mod 125) P125,1 
88 (mod 200) 601 ir od 15) 
188 (mod 200) 7137001 RR? 


For Table A.39, we set 
po = 236661696642275153056980146191674776616380367693641, 


pos = 24106981477091678423113880081946849059226586740161, 
p33 = 263768160996144192120004532942855021486 7605295594581 16494001319, 
p35 = 3822099625066 14718527774493539296028007087595862045098246402975616— 
281748634348102379. 


For Tables A.42 and A.44, we set 
Pi2 = 261116663697161542351918133573442849307, 
P15 = 3933848803218900604400041, 


Piz9 = 3889436310686727916228493162492361601. 
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Table A.39 First part of covering for A + 12-31* 


row congruence prime p; 

1 k =3 (mod 6) c 

2 k =4 (mod 6) 19 

3 k =1 (mod 4) 13 

4 k =8 (mod 16) 17 

5 k =0 (mod 11) 23 

6 k =1 (mod 11) 397 

7 k = 2 (mod 11) 617 

8 | k=3 (mod 11) 150332843 

9 | k=4 (mod 22) 757241 

10 | k=15 (mod 22) 1048563011 

Et k = 8 (mod 33) 650141690025315305584300036801 
12 | k=5 (mod 44) 2729 

13 k = 16 (mod 44) 245911396799577828131028569 
14 | k= 27 (mod 88) 89 

15 k = 71 (mod 88) 41440739086 7564627396249 
16 k = 38 (mod 88) | 12236290645201501169749559350025041 
17 | k = 82 (mod 176) P76 

18 | k=170 (mod 352) 309 

19 | k = 346 (mod 352) P352,1 

20 | k=6 (mod 55) 167767051 

21 k =17 (mod 55) pat 

22 | k = 28 (mod 110) 661 

23 | k =83 (mod 110) 2531 

24 | k =39 (mod 110) 11551 

25 | k =94 (mod 110) Pos 

26 | k =50 (mod 165) 1321 

27 k = 105 (mod 165) Pigs 

28 | k = 160 (mod 330) 2995081 

29 k = 325 (mod 330) P330,1 

30 | k=7 (mod 77) 2927 

31 | k= 18 (mod 77) 23503054499 

32 | k= 29 (mod 77) 16169321243923 

33 | k=40 (mod 77) P33 

34 | k=51 (mod 154) 818819 

35 | k = 128 (mod 154) P35 

36 | k =62 (mod 231) 463 

37 | k = 139 (mod 231) 13780537 

38 | k = 216 (mod 231) 816786763717 
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Table A.40 Second part of covering for A+ 12-31* 


TOW congruence prime Di 

39 | k=73 (mod 385) 2311 

40 k = 150 (mod 385) P385,1 

41 k = 227 (mod 385) P385,1 

42 | k = 304 (mod 770) 2868251 

43 k = 689 (mod 770) P770,1 

44 k = 381 (mod 770) P770,2 

45 k = 766 (mod 1540) P1540,2 

46 k = 1536 (mod 1540) P1540,2 

47 | k =20 (mod 99) 199 

48 k = 53 (mod 99) 991 

49 | k=86 (mod 99) 204733 

50 k = 32 (mod 99) 36093579787 

51 k = 65 (mod 99) 294573316951 

52 | k=98 (mod 99) 215792743120601131 
53 k = 19 (mod 99) 3763784187326467459 
54 k = 52 (mod 99) 869535983092745596321 
55 k = 85 (mod 198) 8713 

56 | k= 184 (mod 198) 430057 

57 k = 31 (mod 66) 67 

58 | k =42 (mod 66) 297991 

59 k = 43 (mod 66) 34731987261785578083133 
60 k = 54 (mod 132) 599382278617 

61 k = 120 (mod 132) 169301958609793153 
62 k = 30 (mod 132) 4451983606421686827580205284201 
63 | k =96 (mod 264) 183516169 

64 k = 228 (mod 264) P264,1 


Table A.41 First part of covering for A +18 - 31* 


row congruence prime p; 
1 k =4 (mod 6) i 
2 k =2 (mod 6) 19 
3 k =2 (mod 4) 13 
4 | k=0 (mod 15) 2521 
5 k =5 (mod 15) 327412201 
6 | k=10 (mod 45) 211 
7 | k =25 (mod 45) 63901 
8 | k=40 (mod 45) 106291 
9 k =1 (mod 17) | 751670559138758105956097 
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Table A.42 Second part of covering for A+ 18 -31* 


row congruence prime p; row congruence prime p; 

10 | &=6 (mod 34) 103 35 | k= 60 (mod 306) 307 

11 k = 23 (mod 34) 6841661642646463343047 36 k = 213 (mod 306) 78198146102753533 
12 k = 11 (mod 51) P12 37 k = 77 (mod 306) 19220209997787857 
13 | k= 28 (mod 51) 1961163283 38 | k = 230 (mod 306) Dane 

14 | k = 45 (mod 102) 2796214962413636917873 39 k = 94 (mod 459) 1131567193 

15 k = 96 (mod 102) | 195333779873358973907838097 40 k = 247 (mod 459) P459,1 

16 | k=16 (mod 68) 1399577 41 | k=400 (mod 459) Pas9,2 

17 k = 33 (mod 68) 224499664484159761 42 k =111 (mod 612) 15913 

18 k = 50 (mod 68) 1682325489672499143634073 43 k = 264 (mod 612) 34273 

19 | k = 67 (mod 136) 137 44 k = 417 (mod 612) P612,1 

21 | k=4 (mod 119) 239 46 | k= 128 (mod 1224) 29841121 

22 k = 21 (mod 119) 2857 47 k = 281 (mod 1224) 209530441 

23 k = 38 (mod 119) 78541 48 k = 434 (mod 1224) 330188689 

24 | k=55 (mod 119) Prii9 49 k = 587 (mod 1224) 1133681368703716513 
25 k = 72 (mod 238) 11781795397 50 k = 740 (mod 1224) | 61410533730626180017 
26 | k=191 (mod 238) 2744226674742050863 51 | k = 893 (mod 1224) poaia 

27 | k=89 (mod 238) P238,1 52 | k = 1046 (mod 1224) P1224,2 

28 | k = 208 (mod 238) P238,2 53. | k = 1199 (mod 2448) 637045489 

29 | k = 106 (mod 357) 637603 54 | k = 2423 (mod 2448) P2448,1 

30 | k = 225 (mod 357) 2112727 tats) k = 145 (mod 1071) 39307843 

31 | k = 344 (mod 357) P357,1 56 | k = 298 (mod 1071) Pio71.4 

33 k = 26 (mod 153) 42834601810502407 08 k = 604 (mod 2142) 03001 

34 | k = 48 (mod 153) Pri53 59 | k = 1675 (mod 2142) 8281014336631 
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Table A.43 Third part of covering for A +18 -31* 


row congruence prime p; row congruence prime p; 
60 | k= 757 (mod 2142) | 17562188355029137 84 k = 323 (mod 680) | 3576121 
61 k = 1828 (mod 2142) P2142,1 85 k = 663 (mod 680) P680,1 
62 | k=910 (mod 2142) p2142,2 86 k =71 (mod 765) Tas 
63 = = 1981 (mod 4284) 5714857 87 | k=156 (mod 765) 6121 
64 = 4123 (mod 4284) P4284,1 88 | k=241 (mod 765) | 1074061 
65 - = 1063 (mod 3213) | 7812268926269851 89 k = 326 (mod 765) | 8197741 
66 | k = 2134 (mod 3213) pai 90 | k=411 (mod 765) | press 
67 | k = 3205 (mod 3213) P3213,1 91 k = 496 (mod 765) P765,2 
68 | &=31 (mod 85) 108971 92 | k =581 (mod 1530) | piso. 
69 | k =36 (mod 85) 391206807721 93 | k = 1346 (mod 1530) | pisso2 
70 k =41 (mod 85) 21736504684553261 94 | k = 666 (mod 2295) 4591 
71 | k=46 (mod 85) Pe 95 | k= 1431 (mod 2295) | 9181 
72 k =51 (mod 170) 232122807601 96 | k = 2196 (mod 2295) | 59671 
73 | k= 1836 (mod 170) P170 97 | k=751 (mod 2295) |  poo951 
74 | k =56 (mod 255) 934831 98 | k = 1516 (mod 2295) | pov95.2 
75 | k=141 (mod 255) P255,1 99 | k = 2281 (mod 4590) | 2519911 
76 k = 226 (mod 255) P255,2 100 | k = 4576 (mod 4590) P4590,1 
78 | k= 146 (mod 340) 3061 102 | k=161 (mod 425) P4252 
79 k = 231 (mod 340) 30941 103 | k = 246 (mod 850) 40801 
80 | k = 316 (mod 340) 2964461 104 | k =671 (mod 850) P850,1 
81 | k =68 (mod 340) 5153381 105 | & =331 (mod 850) | psso2 
82 | k=153 (mod 340) P340,1 106 | k = 756 (mod 1700) | 630701 
83 | k = 238 (mod 340) P340,2 107 | k = 1606 (mod 1700) | 865301 
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Table A.44 Fourth part of covering for A+ 18- 


Res 


row congruence prime p; row congruence prime p; 
108 | k = 416 (mod 1275) P1275,1 132 | k =33 (mod 160) 8084410241 
109 | k = 841 (mod 1275) P1275,2 133 | k = 113 (mod 320) 259201 
110 | k = 1266 (mod 2550) 2501 134 | k = 273 (mod 320) 59826881 
111 | k = 2541 (mod 2550) P2550,1 135 k= 63 (mod 160) P160,1 

112 k = 81 (mod 595) 89278561 136 | k = 143 (mod 160) P160,2 

113 | k = 166 (mod 595) 108837401 137 | k = 73 (mod 320) P320,1 

114 | k= 251 (mod 595) 42498121331 138 | k = 233 (mod 320) P320,2 

115 | k= 336 (mod 595) P15 139 | k = 153 (mod 480) 20641 
116 k = 421 (mod 595) P595,1 140 | k = 313 (mod 480) P480,1 

117 | k& =506 (mod 595) P595,2 141 | k = 473 (mod 480) P4802 

118 | & = 591 (mod 1190) 2381 142 k = 3 (mod 45) 337048683633480845467801 
119 | k = 1186 (mod 1190) P1190,1 143 k = 8 (mod 90) 2065411 
120 k = 2 (mod 5) 11 144 | k=53 (mod 90) | 300392264044249601502733598251 
121 k = 13 (mod 16) 17 145 | k =13 (mod 135) 1344742561 
122 k =5 (mod 64) 4801 146 | k =58 (mod 135) P1351 

123 k = 21 (mod 64) Po 147 | k = 103 (mod 135) 7135.2 

124 k = 37 (mod 128) 207 148 | k = 24 (mod 180) 8728381 
125 | k=101 (mod 128) 641 149 | k = 84 (mod 180) 14398921 
126 k = 53 (mod 128) 2689 150 | k = 144 (mod 180) P80 

127 | k=117 (mod 128) 9601 151 | k=29 (mod 30) 880374069121 
128 k =3 (mod 80) 136046551681 152 k = 9 (mod 60) 61 

129 k = 43 (mod 80) P129 153 | k =39 (mod 60) 25621 

130 k = 23 (mod 160) 380641 154 | k=19 (mod 60) 1529401 
131 | k= 103 (mod 160) 1176641 155 | k=49 (mod 60) 304643210761 


For Tables A.45 and A.46, we set 


Por = 1139174673410619194105634164045563117625339372157650642879811073, 


p39 = 62322419393153627851729037464684263699383389269055382039663, 


par = 147882001432537751112306358052999119715341090542830827. 


Table A.45 First part of covering for A + 20 -31* 


row congruence prime p; 

1 k =5 (mod 6) 7 

2 k =4 (mod 16) 17 

3 | k=0 (mod 7) 917087137 

4 | k=1 (mod 14) 11971 

5 | k=8 (mod 14) 71821 

6 k =2 (mod 21) 43 

7 k =9 (mod 21) 6301 

8 | k=16 (mod 21) 2813432694367 
9 k =3 (mod 28) 29 

10 | &=10 (mod 28) 7253 

il | k=17 (mod 28) 13469 

12 | k= 24 (mod 28) 277739477 

13 | k=4 (mod 35) 319061 

14 | k=11 (mod 35) 203633641 

15 k = 18 (mod 35) 9240957640390889951861 
16 | k=25 (mod 70) 71 

17 | & =60 (mod 70) 149269961 

18 k = 32 (mod 70) 60427990638165876546967231 
19 | k = 67 (mod 140) 281 

20 | k = 137 (mod 140) 106261 

21 | k=5 (mod 42) O11 

22 k = 12 (mod 42) 550469850411853 
23 | k= 19 (mod 84) 163598980 

24 | k=6l (mod 84) 1038949965541 
25 k = 26 (mod 84) 4038949965541 
26 | k =68 (mod 168) 337886977 

97 | k = 152 (mod 168) Par 

28 | k = 33 (mod 210) 13454000701 
29 | k=7 (mod 210) 18396393590821 
30 | k=117 (mod 210) 350121327433921 
31 | k = 159 (mod 210) | 4303134368687145997467938682848881 
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Table A.46 Second part of covering for A+ 20-31* 


TOW congruence prime Di 
32 k = 201 (mod 420) P4209 

33 k = 411 (mod 840) 30241 
34 k = 831 (mod 840) 52081 
30 k = 40 (mod 126) 2143 

36 | k =82 (mod 126) 45376431752737 
37 k = 124 (mod 126) 1652484831253806817 
38 k = 6 (mod 49) 6959 

39 k=18 (mod 49) P39 

AO k = 20 (mod 98) 2932755253 
41 k = 69 (mod 98) Par 

42 k = 27 (mod 147) 883 

43 k = 76 (mod 147) P1471 

44 | k = 125 (mod 147) Diaz 

45 k = 34 (mod 196) 197 

46 k = 83 (mod 196) Prig6 

47 | k =132 (mod 392) 29401 
48 k = 328 (mod 392) 946681 
49 | k=181 (mod 392) P392,1 

50 k = 377 (mod 392) 392.2 

51 | k=41 (mod 245) 491 

52 k = 90 (mod 245) 59848731 
03 k = 139 (mod 245) P245,1 

54 k= 188 (mod 245) P245,2 

55 | k= 237 (mod 490) 1471 

56 k = 482 (mod 490) 32341 
57 k = 48 (mod 294) 159937 
58 k = 97 (mod 294) 9561579721 
59 k = 146 (mod 294) | 64636178950385134849 
60 k = 195 (mod 588) P588,1 

61 k = 489 (mod 588) P588,2 

62 | k= 244 (mod 882) 154351 
63 k = 538 (mod 882) Psg2,1 

64 | k = 832 (mod 882) ee 

65 | k = 293 (mod 1470) 282506071 
66 k = 587 (mod 1470) P1470,1 
67 k= 881 (mod 1470) P1470,2 
68 | k = 1175 (mod 2940) 126421 
69 | k = 2645 (mod 2940) 135241 
70 | k = 1469 (mod 2940) 840841 
71 | k = 2939 (mod 2940) 33183781 
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APPENDIX B 


SAGEMATH CODE 


In this appendix we provide various SageMath programs and functions that were used 
to find primes where a given base has a given order, as well as to check whether primes 
were digitally delicate. We have broken the code into sections to provide additional 


structure. 


B.1 FINDING PRIME DIVISORS OF ®,,(a) 


# Initialize polynomial ring R, order n, and base a. 
# Calculate nth cyclotomic polynomial at a. 

R = ZZ|’x’| 

n=18 # Or the order you’re looking for. 

n 

a=3 # Or the base you’re covering. 


f=R. cyclotomic_polynomial(n) (a) 


# Check if largest prime divisor p of n divides f. 
# If so, replace f by f/p. 
# In each case, check if f is a prime power. 
fp i tee Ee 
p=max(n. prime_ factors ()) 
if f%p!=0: 


f.is_prime_power () 
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else: 
print ("p divides .") 
f=Integer(f/p) 


f.is_prime_ power () 


# Search for small primes congruent to 1 mod n that divide f. 
# Terminate if all primes found (primeprod=f). 
# To look for more primes, change 10°7 to something larger. 
primes=|| 
primeprod=1 
KS 
while primeprod<f and k<min(107~7,(f—1)/n+l): 
if £%(n*k+1)==0 and (nxk+1).is_prime(): 
print (n*k+l1) 
primes. append (n*k+1) 
primeprod=primeprod *(nx*xk+1) 
k=k+1 


# Divide f by primes found above. 
# Check if quotient is a prime power. 
for prime in primes: 
f=Integer(f/prime) 
while gcd(f, prime) >1: 
f=Integer(f/prime) 
print (prime) 
Lf taal 


print ("No primes left .") 
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else: 
print (f) 


f.is_prime_power () 


B.2 CHECKING WHETHER PRIMES ARE DIGITALLY DELICATE 


def baserep(n,b): 
# Returns representation of integer in a particular base. 
# Input: Integers n, b. 
# Output: List representing n in base b. 
# Base case: n less than b. 
if n<b: 
return [n] 
# Recursive step: Peel off units digit , truncate n. 


else: 


return [n%b] + baserep ((n—n%b) /b,b) 


def digdelcheck(p, b): 
# Checks if an integer is digitally delicate in a given 
base. 
# Input: Integers p, b. 
# Output: Boolean boolvar. 
# Initialize return variable , express p in base b. 
boolvar=True 


pbaseb = baserep(p,b) 


# Initialize loop variables. 


i=0 
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k=0 
# Iterate over all digit changes. 
while i!=len(pbaseb) and k!=b and boolvar: 
# Terminate if a prime is found. 
if (Integer(p + (k — pbaseb|[i])*b”i).is_prime() or p 
+ (k — pbaseb[i])*b7~i==1) and Integer(p + (k — 
pbaseb[i])*b~i)!=p: 
print(p + (k — pbaseb|i])*b7~i) 
boolvar=False 
# Increase digit /power as necessary. 
IE ee 
k=k-+1 
PE kb 
if i<len(pbaseb) : 
k=0 


i=i+l 


# Return conclusion. 


return boolvar 


def widedigdelcheck(p,b, limit =100): 
# Checks if an integer is digitally delicate in a given 
base out to a specified number of leading zeros. 
# Input: Integers p, b, limit. 
# If no value is chosen for limit , program checks 100 
leading zeros. 


# Output: Boolean boolvar. 
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# Initialize return variable , express p in base b. 


boolvar=True 


pbaseb = baserep(p,b) 


# Initialize loop variables. 
ee 
k=0 
# Iterate over all digit changes. 
while i!=len(pbaseb) and k!=b and boolvar: 
# Terminate if a prime is found. 
if (Integer(p + (k — pbaseb|[i])*b”i).is_prime() or p 
+ (k — pbaseb[i])*b~i==1) and p + (k — pbaseb[i]) «x 
be ude= (ie 
boolvar=False 


# Increase digit /power as necessary. 
if k<b: 
k=k+1 
ip k=: 
if i<len(pbaseb): 
k=0 
i=i+l 
# If p is digitally delicate , check leading zeros. 
if boolvar: 
# Initialize loop variables. No longer need to 
consider delta = 0. 
i=len (pbaseb )+1 
ke. 
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# Iterate over all digit changes. 
while i<len(pbaseb)+limit and k!=b and boolvar: 
# Terminate if a prime is found. 
if Integer(p + (k)*b7i).is_prime(): 
boolvar=False 
#print(ptk*b7i, k, i) 
# Increase digit /power as necessary. 
if k<b: 
k=k+1 
if k==b: 
if i<len(pbaseb )+limit : 
kel 
I=) 
# Return conclusion. 


return boolvar 


def leadingzeroscheck(p,b, limit =100): 
# Checks whether leading zeros of an integer can be 
changed without obtaining a prime. 
# Input: Integers p, b, limit. 
# If no value is chosen for limit , program checks 100 
leading zeros. 
# Output: Boolean boolvar. 


# Initialize return variable , express p in base b. 


boolvar=True 


pbaseb = len(baserep(p,b) ) 
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# Initialize loop variables. 
i=pbaseb+1 
reall 
# Iterate over all possible digit changes. 
while i<pbaseb+limit and k!=b and boolvar: 
# Terminate if a prime is found. 
if Integer(p + (k)x*xb-i).is_prime(): 
boolvar=False 
print (ptk*b°i, k, i) 
# Increase digit /power as necessary. 
if k<b: 
k=k+1 
it~ k= 
if i<pbaseb+limit: 
ket 
i=i+l 
# Return conclusion. 


return boolvar 


B.3 CHECKING WHETHER A SYSTEM OF CONGRUENCES IS A COVERING 


# List congruences k = a (mod m) in two lists. 


# moduli is the list of m. 


# kay is the list of congruence classes a. 


modult S135/°64. 9 18) 18h. 4p 12, BO, 8 72, 144ks, 144, 216; 


216, 216, 108, 216, 216] 


kay = [0, 1, 2, 8, 14, 0, 10, 5, 3, 17, 53, 125, 23, 95, 167, 


71, 143, 215] 
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# Calculate least common multiple of moduli. 
M = Iem( moduli) 
M 


# Create list of residues to check. 
tocover = [i for i in range(M)|] 
# For each residue, check if some congruence covers the 
residue. 
for j in range(len(moduli) ): 
for k in set(tocover): 
if kZ*%moduli| j|==kay [j ]: 
tocover .remove(k) 

# Display what is left to cover. 
tocover 
# Display percentage of residues left. 
float (len (tocover ) /M) 
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